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Mirror matter is a stable self-coUisional dark matter candidate. If parity is a con- 
served unbroken symmetry of nature, there could exist a parallel hidden (mirror) sector 
of the Universe composed of particles with the same masses and obeying the same phys- 
ical laws as our (visible) sector, except for the opposite handedness of weak interactions. 
The two sectors interact predominantly via gravity, therefore mirror matter is naturally 
"dark". Here I review the cosmological signatures of mirror dark matter, concerning 
thermodynamics of the early Universe, Big Bang nucleosynthesis, primordial structure 
formation and evolution, cosmic microwave background and large scale structure power 
spectra. Besides gravity, the effects on primordial nucleosynthesis of the kinetic mixing 
between photons and mirror photons are considered. Summarizing the present status 
of research and comparing theoretical results with observations/experiments, it emerges 
that mirror matter is not just a viable, but a promising dark matter candidate. 
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1. Mirror particles 

In 1956, Lee and Yang proposed the non-parity of weak interactions in their fa- 
mous paper entitled Question of Parity Conservation in Weak Interaction^. In the 
same article they mentioned also a possible way-out to restore the parity, writing 
the following: "// such asymmetry is indeed found, the question could still be raised 
whether there could not exist corresponding elementary particles exhibiting opposite 
asymmetry such that in the broader sense there will still be overall right-left sym- 
metry. If this is the case, it .should be pointed out, there must exist two kinds of 
protons PR and PL, the right-ended one and the left-ended one." These sentences 
can be considered the start of the mirror matter adventure. 

Subsequently experiments confirmed that the weak interactions indeed violate 
the parity symmetry, and the interactions of the fundamental particles arc non- 
invariant under mirror reflection of the coordinate system. The weak nuclear inter- 
action is the culprit, with the asymmetry being particularly striking for the weakly 
interacting neutrinos. For example, today we know that neutrinos only spin with 
one orientation. Nobody has ever seen a right-handed neutrino. 
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Nevertheless, ten years later Kobzarev, Okun and Pomeranchull^ discussed var- 
ious phenomenological aspects of the idea suggested by Lee and Yang, a parallel 
hidden (mirror) sector of particles which is an exact duplicate of the observable 
particle sector. 

In 1974 PavsicP showed that non-conservation of parity in /3-decay does not 
necessarily imply the mirror asymmetry of basic laws of nature, since by including 
the internal structure of particles involved in /3-decay we can in principle restore 
reflection invariance. What is asymmetric at the /3-decay are initial conditions, while 
the weak interaction itself is symmetric. 

Many years later, in 1991, Foot, Lew and VolkaJ^ presented the old original 
idea of mirror matter in the modern context of gauge theories. They developed the 
model with exact parity (mirror) symmetry between two particle sectors, that are 
described by the same Lagrangian and coupling constants, and consequently have 
the same microphysics0 From the modern point of view, the physics of ordinary 
world is described by the Standard Model SU{3) x SU{2) x U{1) with gauge fields 
(gluons, photons, W, Z and Higgs bosons) coupled to ordinary quarks 'Land leptons 
I. In the minimal parity-symmetric extension of the Standard Model,!^ the group 
structure is G®G', where G ~ SU (3) ® SU{2) (E) U (1) and the prime (') denotes, as 
usual, the mirror sector. The physics of mirror sector is described by the analogous 
gauge symmetry SU{'i)' x 5[/(2)' x U{1)' with the corresponding gauge bosons 
(mirror gluons, mirror photons, W and Z' bosons) coupled to mirror quarks q' and 
leptons V . In this model the two sectors are described by the same Lagrangians, 
but where ordinary particles have left-handed interactions, mirror particles have 
right-handed interactions. 

Today we know mirror matter as a stable self- interacting dark matter candidate 
that emerges if one, instead of (or in addition to) assuming a symmetry between 
bosons and fermions (supersymmetry) , assumes that nature is parity symmetric. 
The main theoretical motivation for the mirror matter hypothesis is that it consti- 
tutes the simplest way to restore parity symmetry in the physical laws of nature, 
since it is a matter of fact that the weak nuclear force docs not exhibit parity sym- 
metry. Thereby the Universe is divided into two sectors with the same particles and 
interactions, but opposite handedness, that are connected by universal gravity. 

In other words, for each type of particle, such as electron, proton and photon, 
there is a mirror twin with the same mass, so that the ordinary particles favour 
the left hand, the mirror particles favour the right hand. If such particles exist in 
nature, then mirror symmetry would be exactly conserved. In fact one could have 



^ In this review I do not consider the so-called "broken" or "asymmetric" mirror models, where 
the parity symmetry is spontaneously broken and the two particle sectors have different weak 
interaction scales and then different microphysics. For an example, see Ref. [5] 

^ Astrophysical constraints on self interactions of dark matter are present in literature, but the 
reader should keep in mind that they are valid only for homogeneous distributions of dark matter 
particles, and are therefore not directly applicable to the mirror matter case, as well as to all 
non-homogeneous dark matter distributions. 
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the parity symmetry as an exact symmetry of exchange of the ordinary and mirror 
particles. Thus, mirror particles are stable exactly as their ordinary counterparts. 

According to this theory the ordinary and mirror particles have the same masses, 
and the three non-gravitational forces act on ordinary and mirror sectors completely 
separately (and with opposite handedness: where the ordinary particles are left- 
handed, the mirror particles are right-handed). For example, while ordinary photons 
interact with ordinary matter, they do not interact with mirror matter. Similarly, the 
"mirror image" of this statement must also hold, that is, the mirror photon interacts 
with mirror matter but does not interact with ordinary matter. The upshot is that 
we cannot see mirror photons because we are made of ordinary matter. The mirror 
photons would simply pass right through us without interacting at all! 

The mirror symmetry requires that the mirror photons interact with mirror 
electrons and mirror protons in exactly the same way in which ordinary photons 
interact with ordinary electrons and ordinary protons. A direct consequence of this 
is that a mirror atom made of mirror electrons and a mirror nucleus, composed 
of mirror protons and mirror neutrons, can exist. In fact, mirror matter made of 
mirror atoms would also exist with exactly the same internal properties as ordinary 
matter, but would be completely invisible to us! 

Besides gravity, mirror matter could interact with ordinary mat ter via some 
other messengers, as the so-called kinetic mixing of gauge bosons,!^^^ unknown 
fields that car ry both ordinary and mirror charges, neutrino-mirror neutrino mass 
mixing,IIIH^ and Higgs boson-mirror Higgs boson mixing jf g^^ch interactions 
exist they must be weak and are, therefore, negligible for many cosmological pro- 
cesses, but could be determinant for the detection of mirror dark matter in some 
contexts. Since photons do not interact with mirror baryons, or interact only very 
weakly, the presence of mirror matter is felt mainly by its gravitational effects, 
which is exactly the definition of "dark matter" I This is a key aspect of the mirror 
scenario, since it predicts the existence of dark matter in the Universe in a natural 
way. 

After the first pioneering works on mirror matter during 80'J^1~ISI22HS1]^ its many 
consequences for particle physics and astrophysics have been studied during the last 
decades. The reader can refer to Ref. |36]for a history of the development of mirror 
matter research before 2006. 

Like their ordinary counterparts, mirror baryons can form atoms, molecules and 
astrophysical objects such as planets, stars and globular clusters. Invisible stars 
made of mirror baryonJ^S^SSl are candidates for Massive Astrophysical Compact 
Halo Objects (MACHOs), which have been o bserved as rn icrolensing events. The ac- 
cretion of mirror matter onto celestial 

objectJ233ll32l35] 

an d th e presence of mirror 

matter inside compact stars, for example neutron stars,'^^^^ could have interest- 
ing observable effects. If the dark matter is composed of mirror matter then it is 



"^The presence of the mirror sector would cause an invisible decay of th e Higg s boson, that could 
be detected in next measurements at the Large Hadron Collider (LHC)! ^'^ ! ^^^ 
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also possible that some mirror matter exists in our solar system and in other solar 
systems. In particular, it can h ave implications for the s tudy of anomalous events 
within our solar system,l3SHS3] dose-in extrasolar planets, ^S^HSU Pioneer spacecraft 
anomalieslSQIpurthermore, the mirror matter hypothesis has been invoked in various 



physical and astrophysical quest ions P^^ ^ The consequences of mirror matter on 
Big Bang nucleosynthesis (BBN) IS5HMI primordial structure formation, cosmic mi- 
crowave background (CMB) and large scale structure (LSS) of the Universe,!^^!"!^ 
have been carefully investigated. For reviews on mirror matter theory and conse- 
quences, see Refs. [Ml [HI USHHl 

All these studies provide stringent bounds on the mirror sector and prove that 
it is a viable candidate for dark matter. In addition, mirror matter provides one 
of the few potential explanations for the recent DAM A/LIBRA annual modulation 
signal ,1^2™^ and possibly other experiments, as suggested for low energy electron 
recoil data from the CDMS collaboration,!^ the two events seen in the CDMSII/Ge 
experiment ,^221 the rising low energy spectrum observed by CoGeNT.^SS 

If there exists a small kinetic mixing between ordinary and mirror photonsj ^ ^ ^^ ^ 
the mirror particles become sort of "millichargcd" particles for the ordinary ob- 
server, which suggests a very appealing possibility of their detection in the ex- 
periments for the direct search of dark matter, as mentioned above. It is also 
extremely interesting that such a kinetic mixing can be independently tested in 
laboratory "table-top" exp eriments for searching the ortliopositronium oscillation 
into its mirror counterpart .I^^^^HSSl The mixing between the ordinary and mirror 
neutrinos coul d p rovide a possible mechanism for the generation of ultra high en- 
ergy neutrinos^ and for t he gamma ray bursts as a result of explosion of the 
mirror supernovae ! ^ O^ hDiMI Explosions of the mirror supernovae and possible en- 
ergy transfer of ordinary supernovae in the mirror sector could provide a necessary 
energy budget for heating the gaseous part of the mirror matter in the galaxies, 
and hence to prevent its collapse to a disk. In this case the mirror matter could 
form spheroidal halos in accord to observations .^1^ In addition, the efficiency of 
mirror supernovae explosions can be indirectly tested in the future detectors for the 
gravitational waves. 

The mirror dark matter paradigm implies a spectrum of dark matter particles 
with known masses, given by the masses of the stable nuclei (H', He', O', 
Fe', According to this theory, mirror gas and dust should constitute the dark 

spheroidal galactic halos of galaxies .^^^^ Inside the galaxy, mirror baryons should 
form astronomical objects and distribute in an complex way, eventually mixed wit h 
ordinary baryons,!^! "^^ and even accumulate into visible compact stars.^^ 



As usual, mirror quantities are denoted with a prime ('). 
" The distribution of these elements in the dark component of the galaxy could be determined, in 
principle, once we know the chemical evolution of mirror nuclei, starting from their genesis in the 
mirror primordial nucleosynthesis and going on through successive stages of stellar formation and 
evolution. At present we still lack the stellar formation, that constitutes a key ingredient in this 
study. 
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It is worthwhile to note that the presence of the mirror sector does not in- 
troduce any new parameters in particle physics (if we neglect the possible weak 
non-gravitational interactions between visible and hidden sectors). Even more im- 
portant, the mirror particles can exist without violating any known experiment or 
observation. Thus, the correct statement is that the experiments have only shown 
that the interactions of the known particles are not mirror symmetric, they have 
not demonstrated that mirror symmetry is broken in nature. 

Given its consistency with experiments and observations, and the unfruitful 
attempts to prove the existence of the other dark matter candidates, scientific com- 
munity is facing an emergent question: "is mirror matter the dark matter of the Uni- 
verse (or at least a significant part of it)?" One possibility to answer this question is 
to look at the cosmological signatures of mirror particles, that I try to summarize in 
this review. Here I consider the cosmological effects of the gravitational interaction 
between ordinary and mirror matter, and I add the effects of photon-mirror photon 
kinetic mixing only when describing the primordial nucleosynthesis. 

Photon— mirror photon kinetic mixing. Among the possible non-gravitational 
interactions between ordinary and mirror particles, photon-mirror photon kinetic 
mixing plays an important role, since it is related to the interpretation of dark 
matter direct detection experiments. _______ 

In fact, it has been shown in Ref. [551 up-dating earlier studies, ^*^ ! ^^!^^ ! that the 
mirror dark matter candidate is capable of expl aining th e positive dark matter sig- 
nal obtained in the DAMA/LIBRA experiment j-^^^^'^l-^^ SBl while also being consistent 
with the results of the other direct detection experiments.l^SHSSlThe simplest possi- 
bility sees the mirror particles coupling to the ordinary particles via renormalizable 
photon-mirror photon kinetic mixing^! (such mixing can also be induced radiatively 
if heavy particles exist charged under both ordinary and mirror U{l)en^ with 
Lagrangian 

r - 1 Ff^" F' fll 

'-mix — 2 ^ /Ji^ ' y'-j 

where F^''' ^ df'A'' - d^Af" and F'f"" = df'A"' - d'^A''' are the field strength tensors 
for ordinary and mirror electromagnetisms, and e is a free parameter. This mixing 
enables mirror (ordinary) charged particles to couple to ordinary (mirror) photons 
with charge eqe, where q = —1 for e' , q = +1 for p' , etc. 

Thus, the effect of this mixing is equivalent to that of very tiny electric charges, 
the so-called "millicharges" In practice, if mirror and ordinary photons can ki- 
netically mix, the mirror (ordinary) charged particles would behave as millicharges 
in the ordinary (mirror) sector. The debate on the existence of fractional electric 
charges in the Universe is still open, with upper bounds that becomes more and 
more stringent. These bounds, coming from laboratory or astrophysics, are corn- 



Given the current upper bounds, the correct name should be "nanocharges" , but the term 
"millicharges" is still commonly used for historical reasons. 
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putcd for generic millicharges. Therefore, one should be eareful to apply them to 
mirror particles, since they can depend on the details of the model, for example the 
spatial distribution of particles or other interactions concomitant with the kinetic 
mixing of photons. 

An important experimental consequence of Eq. ([l} is the mixing of orthopositro- 
nium with mirror orthopositronium, leading to oscillations between these states in 
a vacuum experiment. The subsequent decays of the mirror state are invisible, re- 
snhing in an effective increase of the decay rlte.S An eventual discrepancy between 
the theoretical predictions and the experimental measurements ma y in fact be re- 
solved by this mirror world mechanism! ^'^ ! ^^ A new measurenientff^ performed 
in 2003 solved the previous longstanding discrepancy, so that currently there is 
no e vid ence of the mirror world from orthopositronium studies, placing an upper 



limitPSe < 5 x 10 More recently, improved measurements of the invisible decay 
of orthopo sitr oniunj^ lowered the upper limit to e < 1.55 x 10~^. A new proposed 



experimeniPfil promises the expected sensitivity in mixing strength to be e ^ 10"'^, 
which is in a region of parameter space of great theoretical and phenomenologi- 
cal interest. In fact, the mirror matter interpretation of the recent DAMA/LIBRA 



observations of the annual modulation signaP"' requires just e ^ 10 . wh ich is 
currently consistent with laboratory and astrophysical constraints! ' '^U'J° I 1U'^ I 



2. Mirror baryons as cosmological dark matter 

In the previous section we have seen that, if we neglect possible but very small 
non-gravitational interactions, the presence of the mirror sector does not introduce 
any new parameters in particle physics. 

However, the fact that microphysics is the same does not mean that also macro- 
scopic realizations of the two particle sectors should be the same. The different 
macrophysics is usually parametrized in terms of two "cosmological" free parame- 
ters: the ratio x of temperatures of two sectors, in terms of temperatures of the or- 
dinary and mirror photons in the cosmic background radiation; the relative amount 
/3 of mirror baryons compared to the ordinary ones. Indeed these two parameters 
are not completely free, since there are astrophysical bounds, that become more 
and more stringent with the progress of research, as we will see in the following. 

All the differences with respect to the ordinary world can be described in terms 
of only two free parameters 

/C/\l/3 rp, Ql 

..(I) and P) 

where T (T'), ilf, (f^J,), and S {S') are respectively the ordinary (mirror) photon tem- 
perature, cosmological baryon density (normalized, as usual, to the critical density 
of the Universe), and entropy per comoving volume. These parameters are defined 
in order to be invariant along the evolution of the Universe, if entropies per comov- 
ing volume and baryon numbers are separately conserved in two sectors. As will be 
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more clear in next section, the parameter x is not exactly the ratio of temperatures 
during all cosmic history, but it is a good approximation during most of the time 
for the usually considered values of x. 

Since it interacts only gravitationally with our ordinary sector, mirror matter 
is a natural dark matter candidate. At present there are observational evidences 
that dark matter exists and it's density is about 5 times that of ordinary baryons. 
This is not a problem for mirror scenario, because the same microphysics does not 
imply the same initial conditions in both ordinary an d mirror sect ors. In particular, 

various baryogenesis scenarios with a hidden sectoJSMnnHm 

support a mirror 

baryonic density at least equal to the ordinary oneH The bounds on the parameter 
/3 come from the requirements that mirror baryon density cannot be more than 
that inferred for dark matter, and should be at least of the same magnitude of 
the ordinary baryon density, in order to give a significant contribution to the total 
matter density. This translates into 

1 < < « 5. (3) 

If the mirror sector exists, then the Universe along with the ordinary electrons, 
nucleons, neutrinos and photons, should also contain their mirror partners. One 
could naively think that due to mirror parity the ordinary and mirror particles 
should have the same cosmological abundances and hence the ordinary and mir- 
ror sectors should have the same cosmological evolution. In this case the mirror 
photons 7', electrons e^' and neutrinos v'^ ^ ^ would give a contribution to the 
Hubble expansion rate equivalent to that of the ordinary ones, i.e. in terms of en- 
ergetic degrees of freedom g' = g = 10.75. Since each neutrino family contributes 
Ag = (7/8) • 2, t he m irror contribution is equivalent to an effective number of extra 
neutrino familie Ai V^ = Ni, — 3 = (10.75/1.75) ~ 6.14. However, this is forbid- 
den by current estimateffflSl = 3.68to;7o or = 3.80^°'^°, depending on the 
value used for the neutron lifetime t„. The only known possibility is then to have 



s As far as the two sectors have the same particle physics, it is natural to think that the mirror 
baryon number density nj^ is determined by the baryogenesis mechanism which is similar to the 
one which fixes the ordinary baryon density n;,. Thus, since the mass of mirror baryons is the same 
as the ordinary ones, one could question whether the ratio could be naturally of order 1 or some- 
what bigger. There are several baryogenesis mechanisms as are GUT baryogenesis, leptogenesis, 
electroweak baryogenesis, etc. At present it is not possible to say definitely which of the known 
mechanisms is responsible for the observed baryon asymmetry in the ordinary world. However, it 
is most likely that the baryon asymmetry in the mirror world is produced by the same mechanism 
and moreover, the properties of the B and CP violation processes arc paramctrically the same in 
both cases. 

But the mirror sector has a lower temperature than the ordinary one, and so at epochs relevant 
for baryogenesis the out-of-equilibrium conditions should be easier fulfilled for the mirror sector. 
In particular, we know that in certain baryogenesis scenarios the mirror wor ld gets a larger baryon 
asymmetry than the ordinary sector, and it is pretty plausible that /3 > l.lfiSlxhis situation emerges 
in a particularly natural way in the leptogenesis scenario due to the lepton number exchange from 
one sector to the other, which leads to > rif,, and can thus explain the near coincidence of 
visible and dark components in a rather natural way.[IinHlSl 
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a lower temperature for the mirror sector, so that the contribution of relativistic 
degrees of freedom at BBN epoch is suppressed by a (TVT )^ factor. In this case 
the contribution of the mirror sector translates intJ^SESl AA''^, w 6.14a:;^, where 
X ~ T' /T (for more details see next section). Thus, the conservative upper bound 
AiV^ < 1 implies the upper limit 

x < 0.64(A7V^)i/4 ~ 0.64, (4) 

with rather mild dependence on ATV^, e.g. AN^, < 0.5 implies roughly x < 0.5 
while ANi, < 0.2 implies x < 0.4. Considering the less conservative upper bound of 
Ref. lll6l even values x < 0.7 are acceptable. It is remarkable that the aforementioned 
indirectly measure d va lues of N^, are only marginally consistent with the Standard 
Model expectation) ^^^ ^ and the presence of mirror particles with a lower temperature 
than ordinary ones would improve the consistency. 

Hence, in the early Universe the mirror system should have a somewhat lower 
temperature than ordinary particles. Therefore, even if their microphysics is the 
same, initial conditions are not the same in both sectors, and in the early Universe 
the mirror particles are colder than the ordinary ones. This situation is plausible if 
the following conditions are satisfied. 

(A) At the Big Bang the two systems are born with different temperatures. Ac- 
cording to the inflationary paradigm, at the post-inflationary epoch the mirror 
sector is (re)heated at lower temperatu re than in the or dinary one, which can 
be naturally achieved in certain models. ^^ ^ ^'^ l ^^ l '*^'*^^ ! '*^'*^'^ ! 

(B) At temperatures below the reheating temperature the two systems interact very 
weakly, so that they do not come into thermal equilibrium with each other after 
reheating. This condition is automatically fulfilled if the two worlds communi- 
cate only via gravity. If there are some other effective couplings between the 
ordinary and mirror particles, they have to be properly suppressed. 

(C) Both systems expand adiabatically, without significant entropy production. If 
the two sectors have different reheating temperatures, during the expansion of 
the Universe they evolve independently, their temperatures remain different at 
later stages, T' < T, and the presence of the mirror sector would not affect 
primordial nucleosynthesis in the ordinary world. 

Therefore, the BBN bounds require that at the nucleosynthesis epoch the tem- 
perature of the mirror sector be smaller than that of the ordinary one, T' < T. 
As far as the mirror world is cooler than the ordinary one, cc < 1, in the mirror 
world all key epochs (as are baryogenesis, nucleosynthesis, recombination, etc.) pro- 
ceed in somewhat different conditions than in ordinary world. Namely, in the sector 
world the relevant processes go out of equilibrium earlier than in the ordinary sec- 
tor, which has many far going implications, in particular on the key cosmological 
epochs. In next sections we will see the consequences for the thermodynamics of the 
early Universe, BBN, structure formation and evolution. 
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In the most general context, the present energy density contains relativistic 
(radiation) component fir, non-relativistic (matter) component fi™ and the vacuum 
energy (cosmological term or dark energy) density f^A- According to the inflationary 
paradigm the Universe should be almost flat, 51o = + + ^^A ~ 1, which agrees 
well with the results on the CMB anisotropy. 

Then, in a so-called Mirror Universe^ fir and ri„i represent the total amount of 
radiation and matter of both ordinary and mirror sectors: fir ~ {^r)o + {^r)M and 

flm = {flm)o + {flm)M- 

In the context of our model, as explained in next section, the relativistic frac- 
tion is represented by the ordinary and mirror photons and neutrinos. Using the 
expression for the ordinary energetic degrees of freedom in a Mirror Universe, 
g{T) = g{T){l + x^) (see Sec. 13.21 for details), and the value of the observable 
radiation energy density {i}r)oh^ — 4.2x10"^ (where h is, as usual, the Hubble 
parameter), it is given by 

fir = 4.2 X 10"^ h^^ (1 -I- x'^) ~ 4.2 X 10"^ h'^ , (5) 

where the contribution of the mirror species, expressed by the additional term x^, 
is very low in view of the BBN constraint x < 0.64. As for the non-relativistic 
component, it contains the ordinary baryon fraction flf, and the mirror baryon 
fraction OJ, = /3J7b, while the other types of dark matter could also be present. 
Obviously, since mirror parity doubles all the ordinary particles, even if they are 
"dark" (i.e., we are not able to detect them now), whatever the form of dark matter 
made by some exotic ordinary particles, there will exist a mirror partner made by 
the mirror counterpart of these particlesi Thus we consider a matter composition 
of the Universe expressed in general by 

n„, = Ob + O't, + flDM = ^bil ^DM , (6) 

where the term floM includes the contributions of any other possible dark matter 
particles but mirror baryons. 

Cosmological key epochs. As shown in Refs. l65 | [70H73l due to the temperature 
difference between the two sectors, the cosmological key epochs take place at dif- 
ferent redshifts, and in particular they happen in the mirror sector before than in 
the ordinary one. 



^ The expression "Mirror Universe" is perhaps misleading, since one could think that there is 
another Universe, while it is just one, but made of two particle sectors, ordinary and mirror. 
Nevertheless, this expression is sometimes used to shortly refer to this scenario. 
' In the context of supersymmetry, the cold dark matter (CDM) component could exist in the form 
of the lightest supersymmetric particle (LSP). It is interesting to remark that the mass fractions 
of the ordinary and mirror LSP are related as f^^^p — x^lsp- In addition, a significant HDM 
component f2^ could be due to neutrinos with order eV mass. The contribution of the mirror 
massive neutrinos scales as n[, = x^Qi, and thus it is irrelevant. In any case, considering the only 
CDM component, which is now the preferred candidate, we can combine both the ordinary and 
mirror components, since their physical effects arc the same. 
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The relevant epochs for the cosmic structure formation are related to the matter- 
radiation equality (MRE) Zoq, the matter-radiation decouplings (MRD) Zdec and 
■^dcc "^^"^ plasma recombinations in both sectors, and the photon-baryon 

equipartitions z^-y and z[^^. The MRE occurs at the redshift 

1 + .., = ^. 2.4. 10*5^, (7) 

which is always smaller than the value obtained for an ordinary Universe, but ap- 
proximates it for low X (see Fig.[T|). If we consider only ordinary and mirror baryons 
and photons, we find 

where the baryon and photon densities refer to the redshift z. This implies that, with 
the addition of a mirror sector, the matter-radiation equality epoch shifts toward 
earlier times as 

The MRD, instead, takes place in every sector only after most electrons and 
protons recombine into neutral hydrogen and the free electron number density Ue 



J In the most general case, where there is also some other dark matter component, the parameter 
/9 could be the sum of two terms, /3 -|- fioM- 



Fig. 1. The mirror photon decoupling redshift 1 -|- z^^^ as a function of x (thick solid). The 
horizontal thin solid line marks the ordinary photon decoupling redshift 1 + z^^^ = 1100. We 
also show the matter-radiation equality redshift 1 -|- Zeq (dash) and the mirror Jeans-horizon mass 
equality redshift 1 + z'^ (dot) for the case 

r2„i/i2 = 0.135 (see Sec. 15. 1.511 . The shaded area x > 0.64 

is excluded by the BBN limits. 
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diminishes, so that the interaction rate of the photons = n^ax = X^rin^aT drops 
below the Hubble expansion rate H{T), where <jt is the Thomson cross section, 
Xf, = n(./ni, is the fractional ionization, and r/ = n}j/n^ is the baryon to photon 
ratio. In the condition of chemical equilibrium, is given by the Saha equation, 
which for ^ 1 reads 

X.«(l-r4)V2^(^)"'^'e-^/2^ (10) 

where B = 13.6 eV is the hydrogen binding energy and Y4 is the "'He abundance. 
Thus we obtain the familiar result that in our Universe the MRD takes place in the 
matter domination period, at the temperature Tdcc — 0.26 eV, which corresponds 
to redshift 1 + Zdcc = Tdcc/7o ~ 1100. 

The MRD temperature in the mirror sector T^^^ can be calculated following the 
same lines as in the ordinary one. Due to the fact that in either case the photon 
decoupling occurs when the exponential factor in Eq. (|10p becomes very small, we 
have T^^^ ~ Tkcc, up to small corrections related to r/, Y( which are respectively 
different from rj, Y4. Hence, considering that T' = x ■ T, wc obtain 

1 + z^,, ~ x-i(l + Zdcc) ^ 1.1 • lO^x-i , (11) 

so that the MRD in the mirror sector occurs earlier than in the ordinary one. 
Moreover, comparing Eqs. ([7]) and (jlip . which have different dependences on x, we 
find that, for x smaller than a typical value Xcq expressed by 

Xeq » 0.046(a„/i')"' , (12) 

the mirror photons would decouple yet during the radiation dominated period (see 
Fig-HJ- This quantity has a key role in structure formation with mirror dark matter. 
In fact we expect that for values x < Xeq the evolution of primordial perturbations 
in the linear regime is practically identical to the standard CDM case. Assuming, 
e.g., the value ilmh'^ = 0.135, we obtain Xcq ~ 0.34, which indicates that below 
about this value the mirror decoupling happens in the radiation dominated period, 
with consequences on structure formation (as we will see in the following sections 
and as shown in Refs. [65l [69H73|) . 

The redshifts of photon-baryon equipartitions for ordinary and mirror particles 
are related by 

l + 4. = gf-^ = (l + ^..)|>l + ^..- (13) 

This means that also photon-baryon equipartition happens in the mirror sector ear- 
lier than in the ordinary one. 

In Sec. Owe will study the important consequences of these effects on structure 
formation. 
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3. Thermodynamics of the early Universe 
3.1. Thermodynamical equilibrium 

We extend the standard theory of thermodynamics in the early Universe in order 
to take the existence of mirror particles into account. 

We consider the Universe as a thermodynamical system composed of different 
species (ordinary and mirror electrons, photons, neutrinos, nuclcons, etc.) which, in 
the early phases, were to a good approximation in thermodynamical equilibrium, 
established through rapid interactions, in the two sectors separately. 

We assume, as usual, that the Universe is h omogene ous , and the chemical po- 
tentials of all particle species A are neglig 

^^g| 118 | 119 | ^ The latter 

assumption implies the conservation of the total entropy of the Universe S. There- 
fore we can use the equilibrium Bose-Einstein or Fermi-Dirac distribution functions 
and calculate the energy density p, the pressure p, and the entropy density s for 
every particle species in thermal equilibrium 

'^^^^-2^L exp(i?/r)±l 

Pa{T) = / , , dE , 15) 

67r^ Jm. ^ exp {E/T) ± 1 

/^N PAiT) + pA{T) 
SA[T) = , (16) 

where gA is the spin-degeneracy factor of the species A, the signs -I- and - correspond 
respectively to fermions and bosons, and E = -y/p^+m^, with p the momentum. 
We can use the usual parametrization of the entropy density 

siT) = ^qiT)T\ (17) 



where 



1{T)^ E E 9f{T)(^y (18) 



T J 8 ^ ' ' V T 

fcrmions 



is the number of effective entropic degrees of freedom (d.o.f.), that group together 
the d.o.f. for all bosons (gb) and fermions (<?/). An analogous formalism can be used 
for the total energy density p, which can be parametrized as 



7r2 



where 



PiT) = ^ 9iT) , (19) 



9iT)^ E + J ^ 9AT)(^) (20) 

bosons fcrmions 



is the number of effective energetic degrees of freedom. 
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3.2. A useful simple approximation 



Once the ordinary and mirror systems are decoupled already after reheating, at 
later times t they will have different temperatures T{t) and T'{t), and so different 
energy and entropy densities 

Pit) = ^5(T)r4 , p'it) = ^ff'(r')r'4 , (21) 
s{t) = ^q{T)T' , At) = |^g'(r')T'3 . (22) 

The factors g, q and g' , q' account for the effective number of the degrees of freedom 
in two systems, and in general can be different from each other. 

The Hubble expansion rate is determined by the total energy density p = p + p' , 
H = ^ [^i: /'i)Gp. Therefore, at a given time t in a radiation dominated epoch we 
have 

Hit) = 1 . 1.66V^^ = (23) 

in terms of ordinary and mirror temperatures T{t) and T'{t), where 

giT) = g{T){l + ax^) ^ g{T){l + x^) , (24) 



-g'iT') = g'{T') (l + ^) ^ d'iT') (l + x'^) 



(25) 



and A/pi ~ 1.22 x 10^^ McV is the Planck mass. Here the factor a{T,T') = 
[g'{T')/g{T)] ■ [q{T) / q' {T')]'^/^ takes into account that for T' ^ T the relativis- 
tic particle contents of the two worlds can be different. However, except for very 
small values of x, during most of the time in history of the Universe a useful simple 
approximation is a 1. In particular, in the modern Universe we have a(ro, Tq) = 1, 
q{To) ~ q'iTg) ~ 3.91, and x ~ Tq/Tq, where Tq,Tq arc the present temperatures 
of the ordinary and mirror relic photons. 

As far as cc^ ^ 1 , the effective degrees of freedom and the Hubble expansion rate 
are dominated by the density of ordinary particles and the presence of mirror sector 
has negligible effects on the standard cosmology of the early ordinary Universe. The 
opposite happens instead for the mirror sector, that is largely influenced by the 
density of ordinary relativistic particles, as evident from the factor a;~* in Eq. ([25]). 
This makes the cosmology of the early mirror world different from the standard one 
as far as the crucial epochs like baryogcnesis, nucleosynthesis, etc. are concerned. 
Any of these epochs is related to an instant when the rate of the relevant particle 
process r(T), which is generically a function of the temperature, becomes equal to 
the Hubble expansion rate H{T). Obviously, in the mirror sector these events take 
place earlier than in the ordinary sector, and as a rule, the relevant processes in the 
former freeze out at larger temperatures than in the latter. 

It is useful to note that, due to the difference in the initial temperature conditions 
in the two sectors, reactions at the same temperature = occur at different 



January 13, 2013 14:54 WSPC/INSTRUCTION FILE mir-cosmo-rcvicw 



14 Paolo Ciarcelluti 

times t'^_ = a;^i*, which implies different rates of the Hubble expansion (in particular 
H{t'^) > H{t^,), given x < 1), while reactions at the same time t^, = t'^ occur 
at different temperatures = xT,. In particular, this behaviour has important 
consequences on Big Bang nucleosynthesis, that are shown in details in Sec. |31 

3.3. A more accurate study 

During the expansion of the Universe, the two sectors evolve with separately con- 
served entropies, which means that the ratio of entropy densities is also conserved, 
so that in general we can approximate 



This approximation is valid only if the temperatures of the two sectors are not 
too differen t, o r otherwise if we are far enough from crucial epochs, like e^-e~ 
annihilation!^ If, instead, we are studying the range of temperatures interested 
by this phenomenon, the approximation (|26p is no more valid, and a more accu- 
rate study is required for different reasons. First of all, we need a detailed study 
of thcrmodynamical evolution of the early Universe in order to obtain a reliable 
description of the thermal cosmic history in presence of mirror dark matter. Sec- 
ond, an accurate study of the trend of number of degrees of freedom is necessary 
for any future simulation of nucleosynthesis of primordial elements in both sectors, 
that may be compared with observations. In particular, it is even more important 
in view of the recently proposed interpretation of the DAMA/LIBRA experiment in 
terms of interactions with mirror heavy elements.'^ Third, claims for variations of 
the effective number of extra neutrino families, computed at BBN (^ 1 MeV) and 
CMB formation (< 1 eV) epochs, require investigations for possible variations of 
these numbers as consequences of physics beyond the Standard Model (see Ref . 11201 
and references therein for previous works). 

3.3.1. Equations 

We calculate the equations which link the ordinary and mirror sector temperatures 
and thcrmodynamical quantities; then we solve these equations numerically. Once 
the temperatures are known, it is possible to obtain the exact total number of 
d.o.f. in both sectors, which, as common in literature, can be expressed in terms of 
extra-neutrino number. We report in Sec. 13.3.21 the results of these calculations. 

The presence of the other sector leads in both sectors to the same effects of 
having more particles. As already stated, we do not take interactions between the 
two sectors besides gravity into account. This implies that the entropies of the two 
sectors are conserved separately, and the parameter x is constant during the cosmic 
evolution 




(26) 




(27) 
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where a is the scale factor of the Universe. 

As aheady stressed, ordinary and mirror sectors have the same microphysics; 
therefore we can assume, to a first approximation, that the neutrino decouphng 
temperature]^ T^^, is the same in both of them, that is Td^ = T'j^^, . But the tem- 
peratures in the ordinary sector when the ordinary and mirror neutrino decouphngs 
take place are different, Tjji, ^ Tu^/, with Tu^ < Td^> since a; < 1. Therefore, at 
T 3> Tdi, (T'j^^), ordinary (mirror) neutrinos are in thermal equilibrium with the 
ordinary (mirror) plasma and their temperature is = T (T^ = T'), while after 
decoupling it scales as a^^. 

In each sector, shortly after the u decoupling, the e~^-e~ annihilate because 
the temperature becomes lower than 2me, which is the threshold for the reaction 
7 ^ e^e~ in both ordinary and mirror worlds. Thus electrons and positrons transfer 
their entropy to the corresponding photons, which become hotter than neutrinos. 

This fact will be used together with the entropy conservations (total and in 
each sector separately) to find the equations that govern the evolution of the mirror 
photon temperature T' as a function of the ordinary one T . In fact, ordinary and 
mirror neutrino decouplings are key events, together with both e+-e~ annihilation 
processes, for the thermodynamics in this range of temperatures. Once we call Td^' 
the ordinary world temperature when the mirror neutrino decoupling takes place, 
we can split the early Universe evolution for temperatures T < 100 MeV into three 
phases. 



Phase T > Tdv'- Photons, electrons, positrons and neutrinos are in thermal 
equilibrium in each sector separately, that is Ty = T^, = T , T,' = T'^ = T'Q Using 
Eqs. ([T4)) - ((20)) for particles in both sectors we are able to calculate only the d.o.f. 
number in ordinary or mirror worlds; but to work the total d.o.f. number out, we 
need the mirror temperature T" as a function of the ordinary one T or vice versa. 

As we neglect the entropy exchanges between the sectors (valid since there are 
only gravitational interactions between them), we can obtain both these functions 
using Eq. p7)) and imposing x =const. in 



(28) 



+ s'^ + s', ^ 




^97 ^ 


7 


y3 


5g ~\~ 5-y "f" Su 






7 

-g9. 


T3 



where = 6 and 



2, while qe{T) stands for 

" Se(T) 



9e(T) = - 



7 2tt' 
l5 



(29) 



^ The assumption that the neutrino decoupling is an instantaneous process taking place when 
photons have the temperature T = Tov introduces just a small error (< 1%), that is negligible 
for the precision required here. 

' For simplicity we write T^, = , T'^, = T'^ , T^, = T!^ = T' , = T. 
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with Se(T) defined in (HH), and we have used = = = T' , ^ = ^ T. 
We have used also that, since particles and physics are the same in both sectors, 
the equilibrium distribution functions, and thus the spin-degeneracy factors, are the 
same: = q'^, q-y q'^, q„ = ql- 

Equation (|28l) can be solved numerically in order to obtain the function T'{T) 
for every T > T^^i . 

Phase Tuv < T < Tn^'- At T ~ Td^i mirror neutrinos decouple and right 
after mirror electrons and positrons annihilate, transferring their entropy only to 
the mirror photons, and hence raising their temperature. In the mirror sector the 
entropies of the system (e^',7') and of v' are conserved separately. Nevertheless, 
ordinary photons and neutrinos still have the same temperature T. Therefore we 
have two equations. The first one comes from the conservation of entropies in the 
mirror sector, so that their ratio is equal to the asymptotic value computed at high 
temperatures, when the e^'-e~' annihilation process had not begun yet 



and its solution, obtained using Eq. (pO)) . gives T' as a function of T. Eq. (|3T|) is 
the same as Eq. 1^ but with T^' ^T!^ ^T' ^T^. 

Phase T < Tov- At T ~ Toi> ordinary neutrinos decouple and right after ordi- 
nary electrons and positrons annihilate. Now also in the ordinary sector the entropies 
of the system (e^,7) and u arc conserved separately; therefore we need one more 
equation to calculate the ordinary neutrino temperature T^, as a function of the 
ordinary photon one T 




(30) 




-qe{T) + q^ + -q^ 




(32) 




(33) 
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T'3 H 


--a T' 3 




■7 

-qe{T) + q^ 







(34) 



Eq. dSll) is the same as Eq. §1^ but with T^^T-y ^T^. 

Once both ordinary and mirror photon temperatures are known, it is straight- 
forward to calculate the total energy and entropy densities; then, reversing Eqs. (|17p 
and (|19[) . we can obtain the entropic {q) and energetic {g) number of d.o.f. Calcula- 
tions have been made for several different values of x, as reported in the following 
subsection. 

3.3.2. Numerical calculations 

The equations we introduced above have been solved numerically. Since a; is a free 
parameter in our theory, several values have been used for it — from 0.1 to 0.7 
with step 0.1 or less. In the extreme asymptotic cases T ^ Td^i ~ T^y/x or 
T < Tannei =i 1 McV, wc cxpcct q{T) ~ q'[T'{T)]; therefore, using Eq. ([Ml), m 
these limits the ratio of mirror and ordinary photon temperatures should be x, that 
is T'/xT ~ 1. Instead, when Tdu' >T> T^nn e± we expect q'[T'{T)] < q{T) because 
the e^-e~ annihilation takes place before in the mirror world, leading to a decrease 
of q'{T') and a corresponding increase of T' in order to keep the ratio of entropy 
densities constant; thus we expect T' /xT > 1. Later on, when the ordinary electrons 
and positrons annihilate, they make even T increase and thus the ratio T' /xT 
decreases to the asymptotic value 1. These remarks have been verified numerically; 
the ratio T' jxT is plotted in Fig. [2] for different values of x. 



1.4 



1.3 



1.1 





x = 0.1 






x = 0.2 






x = 0.3 






X 04 






x = 0.5 






x = 0.6 






x = 0.7 







4 3 2 1 

T [MeV] 

Fig. 2. The ratio T' /xT for several values of x. The asymptotic values of this ratio are 1 both 
for high and low temperatures, as expected. 
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Mirror e+-e~ annihilation happens at higher ordinary temperatures for lower x 
values, raising before the temperature of mirror photons; this results in a shift of 
the peaks of Fig. [2] towards higher T . In addition, for lower x the difference in T 
between the two annihilation processes is higher, so mirror photons have more time 
to raise their temperatures before the ordinary ones start to do the same. This leads 
to the change of the shape of the curves. 

Number of degrees of freedom. Using Eqs. p^ - ((T5)) to reverse ([T71) and (|19p . 

we can obtain the total number of ordinary entropic [q] and energetic {g) d.o.f. at 
any temperatures T . We can apply the same procedure to calculate the standard 
values qstd and ^std, as well as the mirror ones q' and g' , and the influence of a 
sector on the other one. 

To a first approximation, we expect q (g) to have a cubic (quartic) dependence on 
X when the temperature is not close to the v decoupling and the e'^-e~ annihilation 
phases: q = qstdi^ + x^), q' gstd(l + a;"^), g = .gstd(l + x''), g' = gstdi^ + x''^) (see 
Sec. [33). 

In the continuation we present the results of accurate numerical calculations. 
In Table [T] some values are reported for special temperatures and several values of 
x. As expected, the total d.o.f. numbers are always higher than the standard and 
increase with x. Moreover, the mirror sector values are higher than the ordinary 
ones by a factor of order x~^ (for q) or (for g). 

In Fig.[3jA) Qstd, q, q' are plotted for the intermediate value x = 0.5; panel (B) 
of the same figure shows the corresponding values of g. In the figure q' (g') has been 
scaled by a factor x^ {x'^)', in this way the asymptotic values are the same than the 
ordinary sector ones because at the extremes T'/T = x (see Fig. [2]). We can see 
that ^, 5, q' and g' have similar trends, but, due to T' < T, the d.o.f. number in 
the mirror sector begins to decrease before (at higher T). 

In Fig.|4]the values oi g in ordinary and mirror sectors arc plotted in comparison 
with the standard for several values of x (from 0.1 to 0.7 with step 0.1). As in Fig.jSJ 
the mirror values g' have been multiplied by x'^. The predicted quartic dependence 
of g on X at the extremes is proved correct. In panel (A) we note that ordinary 
d.o.f. with X < 0.3 are practically identical to the standard case. In addition, the 
plot shape does not change with x in the ordinary sector, while it does in the mirror 
one. This sector evolves with temperature T' ~ xT < T; therefore, for lower x the 
number of d.o.f. begins to decrease below the asymptotic value at high T before 
than for high x. The change of the shape of the plots in panel (B) is related to the 
same physical processes responsible for the analogous effect present in Fig. [21 that 
is due to the e^-e^ annihilation in the mirror sector. 

Number of neutrino families. We know that the SM contains three neutrino 
species; the possible existence of a fourth neutrino has been investigated for a long 
time, also using BBN constraints. This is why in literature one can often find bounds 
on the number of d.o.f. in terms of effective extra-neutrino number AiV^ = Ni, — 3. 
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14 




4 3 2 1 

T [MeV] 

Fig. 3. Total entropic (panel A) and energetic (panel B) degrees of freedom computed in ordinary 
and mirror sectors and for the standard. The mirror values have been multiplied by (q) and 
(g) to make them comparable with the ordinary ones, since {T' /xT) ~ 1. 




Fig. 4. Total energetic degrees of freedom in the ordinary (panel A) and mirror (panel B) sectors 
computed for several values of x. In panel (A) the standard is shown for comparison. In panel (B) 
the mirror values have been multiplied by a;^ to make them comparable with the ordinary ones, 
since (T'/xT) ~ 1. 
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Table 1. Standard and non-standard total d.o.f. num- 
bers for several x values in both ordinary and mirror 
sectors. Temperatures are in MeV. 



T 


0.005 


0.1 


0.5 


1 


5 


9std 


3.91 


4.78 


10.0 


10.6 


10.75 


9std 


3.36 


4.30 


10.0 


10.6 


10.75 






X = 


0.1 






T' 


0.0005 


0.0107 


0.0675 


0.124 


0.511 


q 


3.91 


4.78 


10.0 


10.6 


10.75 


1 


oyio 


oyio 


4U ( z 




iUU / U 


9 


3.36 


4.30 


10.0 


10.6 


10.75 


a' 


33629 


32894 


30032 


44430 


98436 






X = 


0.3 






T' 


0.0015 


0.0321 


0.170 


0.315 


1.50 




4.015 


4.91 


10.3 


10.8 


11.0 


g' 


149 


149 


261 


347 


406 


9 


3.39 


4.34 


10.1 


10.6 


10.8 


-g' 


418.5 


409 


750 


1082 


1325 






X = 


0.5 






T' 


0.0025 


0.0533 


0.263 


0.508 


2.50 




4.40 


5.38 


11.3 


11.9 


12.1 


q' 


35.2 


35.5 


77.3 


90.5 


96.5 


9 


3.57 


4.58 


10.7 


11.2 


11.4 


9' 


57.2 


56.6 


138 


168 


182 






X = 


0.7 






T' 


0.0035 


0.0733 


0.357 


0.704 


3.50 




5.25 


6.42 


13.5 


14.2 


14.4 


q' 


15.3 


16.3 


36.9 


40.6 


42.0 


9 


4.17 


5.35 


12.4 


13.1 


13.3 


9' 


17.4 


18.5 


47.8 


53.3 


55.5 



In general, the effective number of neutrinos N^, is found assuming that all particles 
contributing to the Universe energy density, except electrons, positrons and photons, 
are neutrinos; in formula that means 



4 



g{T) = g^XT) + + ^ • 2iV,. • ( ^ ) , (35) 



which implies 



N^^m^f^n^.^rj , (36) 

8 ' ^ 

Ni, has been worked out using Eq. (pS)) together with the results of previous numer- 
ical simulations; some data are reported in Table [21 while plots for several x values 
and temperatures from to 3 MeV are shown in Fig. [5l We stress that the standard 
value = 3 is the same at any temperatures, while a distinctive feature of mirror 
scenario is that the number of neutrinos raises for decreasing temperatures. Anyway, 
this effect is not a problem; on the contrary it may be useful since recent cosmolog- 
ical data fits give indications for a number of neutrinos at recent times higher than 
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Table 2. Effective number of neutrinos in the or- 
dinary sector for some special cases. Temperatures 
are in MeV. 

T x = 0.1 x = 0.3 X = 0.5 X = 0.7 



ordinary sector 

0.005 3.00074 3.05997 3.46270 4.77751 

0.1 3.00074 3.05997 3.46244 4.76829 

0.5 3.00074 3.05997 3.40706 4.52942 

1 3.00071 3.05202 3.39166 4.49133 

5 3.00063 3.04989 3.38430 4.47563 



at BBN epoch. In Rcf. [HOlthe authors found = 5.2^11 using recent CMB and 
LSS data. We can use the aforementioned mirror feature together with the results 
of a previous work on CMB and LSS power spectrap^^ where the author studied 
the dependence of the spectra on the parameters x and N^, for a flat Universe with 
mirror dark matter. As can be easily evinced from Figs. 11, 12, 14, 15 of Ref. [74j an 
increase of the effective number of neutrinos is well mimicked by an increase of the 
parameter x, and the amounts of the respective increases are in accordance with 
what is required to justify the data of CMB and LSS. Thus, considering the sum of 
these effects, i.e. the raise in before and after BBN, and the similarity of CMB 
and LSS spectra of mirror models and standard N^, with the ones obtained without 
mirror sector but with larger N^, the mentioned discrepancy naturally disappears. 
Similarly, the effective number of neutrinos in the mirror sector can be worked out 
as 

■(!)'■ 

8 

Once again these values are higher than the ordinary ones, but now by a factor x~'^; 
they have been computed numerically and some special values are given in Table [3] 
For lower x this number can become very high, inducing relevant consequences on 
the primordial nucleosynthesis in the mirror sector, that is highly dependent on this 
parameter. 



Table 3. Effective number of neutrinos in the 
mirror sector for some special cases. Temperatures 
are in MeV. 

T' x = 0.1 x = 0.3 x = 0.5 x = 0.7 







mirror sector 
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28.86 
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61435 
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101.4 


28.66 


5 


61432 


761.4 


101.3 


28.59 
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x = 0.5 
x = 0.6 
x = 0.7 
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Fig. 5. Effective number of neutrino families in the ordinary sector for several values of x. 



Predictions for special temperatures. It is possible to obtain the asymptotic 
values of Ni, at T ^ and T <C 71,^6* in a simple way starting from the 
standard values 



.9std(r > Tdu) = 10.75 , 

gstd{T < Tanne±) - 3.36 . 

Without the mirror sector, we have, as expected 

10.75-2 - - -4 
N.{T » Td.) = ^ § = 3 



• 2 



NIT^^T \ 3.36 - 2 (n\ 

N^(l < Janne±) = — [^-^ I = 3 



(38) 



(39) 



Instead, when the mirror sector is present, we can use the previously mentioned 
quadratic approximation g = .gstd(l + x*) a-t the special temperatures we are con- 
sidering, and hence 



10.75(1 + ^4) -2- --4 



^ = 3 + 6.14 



^.(r«Ta„ne±; 



M^ii±:^.(HV.3 + 7.40.4 



(40) 
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From the above equations we can see that the increase AiV^ is 
AN, = N,{T « Ta„„e±) - N,{T » To,) = 




(41) 



This leads to higher than the standard in the presence of the mirror world and to 
a further growth of this parameter at low temperatures. If we assume a conservative 
limit on the effective number of neutrino families N,{T ^ Td,) < 4, this implies 
X < 0.64 and N,{T <C Tannei) ^ 4.26. For some interesting values of the parameter 
X we obtain 



iV^(r>Ti,^) -4.5 and Ar^(T < Tan„e± ) - 4.8 , 
x^O.6 iV^(r> Tc^) ~ 3.8 and Ar^(T < Tan„,±) ~ 3.96 , 

x = 0.3 iV^(r> Ti5^) ~ 3.05 and N„{T <^ T^,,^^±) ^ 3.06 . (42) 



These rough estimates are in agreement with the asymptotic numerical values com- 
puted before and after BBN; nevertheless, the previous detailed analysis of the 
evolution of this quantity is crucial for studies of the nucleosynthesis in both sec- 
tors. 

4. Big Bang nucleosynthesis 

The BBN provides one of the most stringent tests for physics beyond the Standard 
Model. Presently it has only one free parameter, the baryon to photon ratio rj = 
rib/n-y, also called baryonic asymmetry. Since microphysics is the same, there should 
be a similar primordial nucleosynthesis also for mirror particles, but, as anticipated 
in the previous section, the different initial conditions for the mirror sector have a 
big influence on this process. For mirror BBN there are two parameters: the mirror 
baryon to photon ratio ?/ = n'^^/n'^ and the ratio of temperatures x, that enters via 
the degrees of freedom, which determines the approximate Hubble expansion rate 
at a given temperature. Using Eqs. ([23]) and (f25|) we obtain 



Since n'^ = x^n-y, using the definition /3 = ^[^/^Ib = n'^^/ni, (the second equality is 
given to the fact that the mass of mirror baryons is the same as the ordinary ones) , 
we obtain 77' = f3x~^ii. With our assumptions ^ and ^ on x and /?, the inequality 
77' > ?7 is always verified. 

4.1. A simple analytical model for He' abundance 

It is useful first to recall the standard BBN. The relevant time scales are: the 
freeze out temperature of weak interactions, Tw ~ 0.8 MeV, and the "deu- 
terium bottleneck" temperature, Tn ^ 0.07 MeV {tN ^ 200 s). For T above 



H{T') = 1.66^g'{T'){l + x-^) 



Mpi 



(43) 
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Tw the weak interactions are rapid enough to keep the neutron-proton ratio in 
chemical equilibrium. The neutron abundance l™l X„ — rin/nb, is then given by 
Xn{T) = [1 + exp(Am/T)]^^, where Am ~ 1.29 MeV is the neutron-proton mass 
difference. For T < Tw the weak reaction rate Tw — G\T^ drops below the 
Hubble expansion rate H{T) ~ 5.5 T^/Afpi, the neutron abundance freezes out 
at the equilibrium value Xn{T\Y) and it then evolves only due to the neutron de- 
cay: Xn{t) = Xn{Tw) expi—t/r), where r = 886.7 s is the neutron lifetime. At 
temperatures T > Tn, the process p + n^d + j is faster than the expansion of 
the Universe, and free nucleons and deuterium are in chemical equilibrium. The 
light element nucleosynthesis essentially begins when the system cools down to the 
temperature 

Tn 7 ^ ^ 0.07 MeV , (45) 

-ln(,) + 1.51n(^) 

where Bd = 2.22 MeV is the deuterium binding energy, and niN is the nucleon 
mass. Thus, assuming that all neutrons end up in ^He, the primordial ^He mass 
fraction is 

^ - = = Ta;^ - 0.25 . (46) 

As discussed, the presence of the mirror sector with a temperature T' < T has 
almost no impact on the standard BBN in the limit x < 0.64, which in fact has 
been set by uncertainties of the present observational situation. In the mirror sector 
nucleosynthesis proceeds along the same lines. However, the impact of the ordinary 
world for the mirror BBN is dramatic! 

Therefore, comparing the Hubble expansion rate H{T') in Eq. with the 
weak reaction rate Tw{T') ~ G^T'^ , we find a freeze-out temperature 

T{y^{l + x~'Y/'Tw , (47) 
which is larger than T^/, whereas the time scales as 

(obtained using Eq. (j43p and the relation t cx H^^). In addition, r/ is different from 
ry ~ 5 X 10^^*^. However, since Tn depends on baryon density only logarithmically 



™ The mass fraction or abundance of a certain nuclear species is defined as 

Xa^^ , J2X, = 1 (44) 

where A is the atomic number, the number density, and rij, the baryonic number density of 
the Universe. Usually the abundance of ''He is alternatively refered as y = Y4 = Xjje, so that for 
*He' we may use Y' = = Xho'- 
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(see Eq. (US])), the temperature Tj^ remains essentially the same as T/y, while the 



time t'j^ scales as 



t'r 



N 



Thus, for the mirror ^He mass fraction we obtain 



(49) 



Y' ^ 2X;(4) = 



2 exp 










1 + exp 


Am 


.riv(i + x-4)i/6_ 



(50) 



We see that Y' is an increasing function of x ^ and is always bigger than Y. In 
particular, for a: — > one has Y' 10 Thus, we have reached the important 
conclusion that, if the dark matter of the Universe is represented by the baryons 
of the mirror sector, it should contain considerably bigger fraction of primordial 
*He than the ordinary world. In particular, the helium fraction of mirror matter is 
between 20% and 80-90%, depending on the values of x and 77'. 



4.2. An accurate numerical model 

If we want a complete and accurate study of the primordial production of both 
ordinary and mirror elements, we need to numerically solve the equations governing 
the evolution of nuclides and nuclei. This study is required, since the primordial 
chemical composition is an initial condition of the following evolution of the Uni- 
verse. 

As we have seen, the presence of the mirror sector can be parametrized in terms 
of extra d.o.f. number or extra neutrino families; therefore, since the physical pro- 
cesses involved in BBN are the same in both sectors, it is possible to use and modify 
a standard numerical cod e to work out the light elements production. The choice is 
the Kawano code for BBl4I2II122] gjj^gg [f^ jg ^ well-tested and fast program and its 
accuracy is large enough for our purposes. For the neutron lifetime we consider the 
value T = 885.7 s. 

The number of d.o.f. enters the program in terms of neutrino species number; 
this quantity is a free parameter, but instead of using the same number during the 
whole BBN process, we use as input the variable N^(T,x) numerically computed 
following the procedure described in the previous section. The only parameter of 



° In reality, Eq. JSDJ is not valid for small x, since in this case deuterium production through 
the reaction n + p <-¥ d + j can become ineffective. It is possible to calculate that for x < 0.3 ■ 
(rj' ■ lO^")^^'''^, the rate at which neutrons are captured to form the deuterium nuclei becomes 
smaller than the Hubble rate for temperatures T' > T'^. In this case mirror nucleosynthesis is 
inhibited, because the neutron capture processes become ineffe ctiv e before deuterium abundance 
grows enough to initiate the synthesis of the heavier elements.ES] We have to remark, however, 
that with our assumption /3 > 1, the condition x < 0.3 ■ (??' ■ lO^")"^/'^ is never fulfilled and the 
behaviour of is well approximated by Eq. |I50|I . 
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the mirror sector which affects ordinary BBN is x; the baryonic ratio j3 does not 
induce any changes on the production of ordinary nuchdes, but it plays a crucial 
role for the mirror nuclides production. 

In Table |4] I report the final abundances (mass fractions) of the elements pro- 
duced in the ordinary sector at the end of BBN process (at T '--^ 8 • 10~^ MeV) for a 
final baryon to photon ratio rj = 6.14 • 10^^" and for several x values, and compare 
them with a standard scenario. We can easily infer that for x < 0.3 the light element 
abundances do not change more than a few percent, and the difference between the 
standard and x = 0.1 is of order lO^"' or less. For x > 0.3 the differences become 
more important. 



Table 4. Elements produced in the ordinary sector. The last row includes 
all elements with atomic number larger than 7. 





standard 


X = 0.1 


X = 0.3 


X = 0.5 


X = 0.7 


n,/H 


1.161 


1.161 


1.159 


1.505 


2.044 


P 


0.7518 


0.7518 


0.7511 


0.7463 


0.7326 


D/H (10-5) 


2.554 


2.555 


2.575 


2.709 


3.144 


T/H (10-**) 


8.064 


8.065 


8.132 


8.588 


10.07 


^He/H (10-5) 


1.038 


1.038 


1.041 


1.058 


1.113 


4He 


0.2483 


0.2483 


0.2491 


0.2538 


0.2675 


6Li/H(10-") 


1.111 


1.111 


1.124 


1.210 


1.499 


■'Li/H (10-10) 


4.549 


4.548 


4.523 


4.356 


3.871 


''Be/H (10-10) 


4.266 


4.266 


4.238 


4.051 


3.502 


8Li + /H (10-15) 


1.242 


1.242 


1.251 


1.306 


1.464 



Between the different possible non-standard BBN scenarios, it would be of rel- 
evance a more careful inspection at the effects of the presence of mirror particles 
on the predicted abundances of light elements. A more extended investigation in 
this sense is required in order to evaluate if it could help solving the still present 
problems on standard BBN due to the lithium anomalies. 

Even mirror baryons undergo nucleosynthesis via the same physical processes 
than the ordinary ones, thus we can use the same numerical code also for the mirror 
nucleosynthesis. Mirror BBN is affected also by the second mirror parameter, that 
is the mirror baryon density (introduced in terms of the ratio /? = Cl'i^/Qi, ^ 1 -f- 5), 
which raises the baryon to photon ratio rj' ~ l3x~^7]. 

The results are reported in Table El which is the analogous of Table |4] but for a 
mirror sector with /3 = 5 and (3 — 1. We can see that BBN in the mirror sector is 
much more different from the standard than BBN in the ordinary sector. This is a 
consequence of the high ordinary contribution to the number of total mirror d.o.f., 
which scales as ^ x~^ (while in the ordinary sector the mirror contribution is almost 
insignificant, since it scales as ^ x^). As we expect, at higher x the mirror abun- 
dances are closer to the ordinary ones, because the mirror sector becomes hotter 
and thus more similar to the ordinary one (since T' ^ xT). Moreover, we note that 
in general lower /? implies final mass fractions closer to the standard, again as we 



January 13, 2013 14:54 WSPC/INSTRUCTION FILE mir-cosmo-rcvicw 



Cosmology with Mirror Dark Matter 27 

may expect, since for /3 = 1 the two sectors have the same baryonic densities. Hence, 
numerical simulations confirm the previous simple analytical model; the mirror he- 
lium abundance should be much larger than that of the ordinary helium, and for 
cc < 0.5 the mirror helium gives a dominant mass fraction of the dark matter of the 
Universe. In addition, mirror BBN produces much larger abundances of so-called 
metals (elements heavier than lithium), that have a large influence on the opacity 
of mirror matter, which has an important role in several astrophysical processes, as 
for example mirror star formation. The interesting feature that mirror sector can 
be a helium dominated world has important consequences on mirror star formation 
and evolution,!^ and other related astrophysical aspects. In particular, I recall the 
reader's attention to the fact that the predicted dark matter composition, domi- 
nated by mirror helium, in combination with results on mirror stellar evolution,!^ 
support the proposed interpretatio n of the DA MA/LIBRA annual modulation sig- 
nal in terms of mirror dark matter .I^SlinilinSl 



Table 5. Elements produced in the mirror sector. 





X = 0.1 (/3 


= 5) 


X = 0.3 (/3 


= 5) 


X = 0.5 (/3 


= 5) 


X = 0.7 (/3 


= 5) 


n/H 


5.762 -10- 


-25 


2.590 -10" 


22 


1.840 -10- 


20 


1.726 -10" 


19 


P 


0.1735 




0.3646 




0.4966 




0.5924 




D/H 


1.003 -10- 


-12 


4.838 10- 


9 


6.587 10- 


8 


3.279 -10- 


7 


T/H 


9.679 -10' 


-21 


1.238 10- 


13 


2.108 10- 


11 


3.722 -10- 


10 


3He/H 


3.282 -10- 


-6 


3.740 10- 


6 


4.172 10- 


6 


4.691 -10- 


6 


"He 


0.8051 




0.6351 




0.5035 




0.4077 




•^Li/H 


7.478 -10" 


-21 


1.309 10- 


17 


1.016 10- 


16 


3.361 -10- 


16 


^Li/H 


1.996 -10- 


-7 


3.720 10- 


8 


1.535 -10- 


8 


7.962 -10- 


9 


''Be/H 


1.996 -IQ- 


-7 


3.720 10- 


8 


1.535 -10" 


8 


7.962 -10" 


9 


8Li + /H 


4.354 -10" 


-9 


5.926 10- 


11 


3.827 10- 


12 


3.949 -10- 


13 




X = 0.1 (/3 


= 1) 


X = 0.3 (/3 


= 1) 


X = 0.5 (/3 


= 1) 


X = 0.7 its 


= 1) 


n/H 


8.888 -10" 


-17 


1.915 -10" 


16 


2.058 10- 


16 


2.076 -10" 


16 


P 


0.1772 




0.3675 




0.5028 




0.6017 




D/H 


1.331 -10- 


-6 


7.094 -10- 


6 


1.352 10- 


5 


2.235 -10" 


5 


T/H 


3.068 -10" 


-9 


2.190 10- 


8 


4.358 10- 


8 


7.328 -10- 


8 


^He/H 


5.228 -10" 


-6 


6.880 10- 


6 


8.232 -10- 


6 


9.719 -10- 


6 


*He 


0.8226 




0.6326 




0.4974 




0.3984 




'^Li/H 


8.638 -10- 


-15 


1.660 10- 


14 


1.790 10- 


14 


1.951 -10- 


14 


'^Li/H 


5.712 -10- 


-8 


8.930 -10- 


9 


2.948 -10- 


9 


1.120 -10- 


9 


'^Be/H 


5.711 -10" 


-8 


8.878 10- 


9 


2.891 10- 


9 


1.064 -10" 


9 


^Li + /H 


2.036 -10' 


-10 


2.514 10- 


12 


1.657 10- 


13 


1.814 -10" 


14 



4.3. BBN with photon mirror photon kinetic mixing 

So far we have considered just gravitational effects of the existence of mirror mat- 
ter. Now we study the implications of photon-mirror photon mixing for the early 
Universe. In particular, we will check that this kinetic mixing is consistent with 
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constraints from ordinary Big Bang nucleosynthesis as well as more stringent con- 
straints from cosmic microwave background and large scale structure considerations. 
We will then estimate, under some simple and plausible assumptions, the primor- 
dial He' mass fraction. Besides studies of mirror star formation and evolution, and 
the formation of dark matter structures at non-linear scales, this quantity is also 
important since the mirror dark matter interpretation of the direct detection exper- 
iments depends on it. In fact, the velocity dispersion of the particles in the mirror 
matter halo depends on the particular particle species and on the abundance of 
mirror helium Y' , and satisfies^"^ 

(51) 

"^^ 2 - -Y' 

4 

where the index i labels the particle type {i — e', H', He', O', Fe'...), v^ot ~ 254 km/s 
is the local rotational velocity for our galaxy, and rup is the proton mass. While the 
DAMA experiments turn out to be relatively insensitive to Y' , other experiments, 
such as those using electron scattering, exhibit a grea ter sensitivity to Y\ and might 
ultimately be able to measure this parameter .^^SEll 



4.3.1. Thermal evolution I: an approximate model 

In the mirror dark matter scenario, it is assumed there is a temperature asymmetry 
(T' < T) between the ordinary and mirror radiation sectors in the early Universe 
due to some physics at early times (for specific models, see Refs. [Mll551[Ml[TT31lll7p . 
This is required in order to explain ordinary BBN, which suggests that T'/T < 0.7. 
In addition, analysei^^l"!^ based on numerical simulations of CMB and LSS suggest 
T'/T < 0.3 (as we will see in next sections). However, if photon-mirror photon 
kinetic mixing exists, it can potentially thermally populate the mirror sector. For 
example, Carlson and 

Glasho\\El 

derived the approximate bound e < 3 x 10 ^ from 
requiring that the mirror sector does not come into thermal equilibrium with the 
ordinary sector before BBN. We expect the kinetic mixing to populate the mirror 
sector, but with T' < T . Assuming an effective initial condition T' <^ T, we can 
estimate the evolution of T'/T in the early Universe as a function of e, and thereby 
check the compatibility of the theory with the BBN and CMB/LSS constraints on 
T'/T. 

Photon-mirror photon kinetic mixing can populate the mirror sector in the early 
Universe, via the process e+e^ — > e+'e^'. The cross section for e+e^ — > e+'e^' is 

^ = , (52) 

3s 

where s ~ {2E)'^ is the Mandelstam s quantity. This implies a production rate for 
e' and a generation of energy density in the mirror sector 

^^^'^ '^'^''■{E) ~ 2n,+ (Tn,- • 3.15 T , (53) 



dt J dt 

' generation 
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where rig/ includes both e^' and e~', and (E) ~ 3.15 T for a Fermi-Dirac distribu- 
tfon. 

The e~^',e^' wiU interact with each other via mirror weak and mirror electro- 
magnetic interactions, populating the ^' ,v'^,v'^,v'^, and thermalizing to a common 
mirror sector temperature T' . The energy density in the mirror sector is then ex- 
pressed by the second of Eqs. (EH, with g' = 43/4 = 10.75 for 1 < T'(MeV) < 100. 
Thus, when g' is constant, we can write 

(54) 

dt iQ^ dt ^ ' 

Using the expression (j53p for the rate of mirror energy density generation, and 
taking the ratio T'/T in order to cancel the time dependence due to the expansion 
of the Universe (redshift), we have 

T J 20 1 



_ , «^e^n,+ n,- ^ A T , (55) 

dt 3.15 TT g' 

with A ~ (1.91/7r^)a^ ~ 3.32 • 10^^, and where in the last equalit y we used the 
expression for densities of relativistic fermions in thermal cquilibriuirfH^ 

3 C(3) „3 ^rc:\ 

ne+ = Tie- = ^ ■ 

In radiation dominated epoch we can express the time a m 

t-OMlg-'/^^, (57) 
where Mpi ~ 1.22 x lO^^ y^^y jg ^j^q pianck mass. For g = 10.75, Eq. dST]) implies 

^ = -— g'/^ Mp 1 ~ -5.46 — . (58) 
dt 0.6 ^ A/pi ^ ' 

Using together Eqs. ([55)1 and ([55|. we obtain 



where B = (A/5.46)A'fpi ~ 0.74 • 10^^ MeV. Finally, integrating Eq. with the 
initial condition T' = at T = T^, we derive the expression for the T'/T ratio as a 
decreasing function of T 

T' /I 1 \ 1/4 

-~0.52.104ei/2 ' 



T ' \T{McY) T,(MeV) 

The last equivalence is obtained considering Ti ^ 1 MeV and defining e_g = e/10^^. 
Clearly in all these computations we have neglected the change of the T'/T ratio 
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due to the e+-e~ annihilation processes in both sectors (as described in detail in 
Sec. [3]), and the transfer from the mirror sector to the ordinary one (since T' <T). 
Note that the T'/T ratio freezes out when T <2 since the number density of 
becomes Boltzmann suppressed and the process e^e~ — > e^' e~' can no longer 
effectively heat the mirror sector. In addition, after e^'-e^' annihilation the effective 
number of mirror degrees of freedom decreases, and thus we estimate 

rpl 

— ~ 0.2 eL(f for T < 1 MeV . (61) 

We plot this expression in Fig. [7] together with the related He' mass fraction (see 
Sec. Constraints from ordinary BBN suggest that SN^ < 0.5 T'/T < 

0.6. A more stringent constraint arises from CMB and LSS, which suggeslP^H-^ 
T'/T < 0.3 and implies, from Eq. dH]), that e < 3 x 10"^. In the limit T'/T ^ 1 
we find the mentioned old upper bound obtained in Ref. j8] However, this is an 
extrapolation beyond the validity of approximations, since, when T' approaches 
T, the transfer of energy from the ordinary to the mirror sector clearly becomes 
more and more inefficient. Nevertheless, far from this limit, this simple analysis 
shows that a photon-mirror photon kinetic mixing of strength e ^ 10~^ (that is 
able to explain the DAMA/LIBRA annual modulation signal consistently with the 
results of the other direct detection experiments) is consistent with constraints from 
ordinary BBN and CMB/LSS data. 



4.3.2. Thermal evolution II: a more accurate model 

A more precise result for the thermal evolution of the mirror sector can be obtained 
using a more accurate model for the transfer of energy from the ordinary to the 
mirror sector, as in Ref. 11081 In this case it is not possible to solve analytically, but 
we need to compute a numerical solution. 

The generation of energy density in the mirror sector due to the process e'^e^ — > 

e IS now 

dp' 

— = n^+ne-{avM0iS) , (62) 

where £ is the energy transferred in the process, VM01 is the M0ller velocity (see e.g. 
Ref.HMl), and again n^- ~ ~ {3C{3)/2n^)T^ . 

It is useful to consider the quantity p' /p, in order to cancel the time depen- 
dence due to the expansion of the Universe. Using expressions (f2T|) and the time- 
temperature relation (|57p with g ~ 10.75, we find that 

P ) _ -n^-n^+{<TVM0\£) 0.6A/pi 
dT ~ TT^gT* ^T3 ■ ^^'^> 

30 
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Let us focus on {(7Vm0i£)- This quantity is 

{avM0i£) 7 J J , (64) 

J i + e^./Ti + ,g./T dVd^P2 

where we have neglected Pauli blocking effects. If one makes the simphfying as- 
sumption of using Maxwelhan statistics instead of Fermi-Dirac statistics then one 
can show (see the Appendix of Rcf. llOSp that in the massless electron limit 

{(TVM01S) = , (65) 



and Eq. (|63|) reduces to 



^ (66) 



dT r2 



^^27C(3)WMp.^ (67) 



where 



This is essentially the same as Eq. (j59p obtained in the previous approximation. Note 
that the e^' will thermalize with 7'. However, because most of the e*' are produced 
in the low T' <b MeV region, mirror weak interactions are too weak to significantly 
populate the v'e^^^r (i-^- one can easily verify a posteriori that the evolution of T'/T 
for the parameter space of interest is such that GpT'^ <C {y/gT^ /O.SMpi). Thus to 
a good approximation the radiation content of the mirror sector consists of e^',7' 
leading to g' = 11/2 and hence p'/p = {g' / g){T'^ /T^), with g'/g « 22/43. 

Assuming the initial condition T' = at T = T^, Eq. has the analytic 
solution 



r fg 



T \g' 



1/4 r 



1/4 



1 1 

[f ' T,_ 

Let us now include the effects of the electron mass. With non-zero electron 
mass, the evolution of T'/T cannot be solved analytically, but Eq. (|63)) can be 
solved numerically. Note that the number density is 



1 + exp{E/T) 
and, as discussed in the Appendix of Ref. 11081 

1 



I ^ " dE (69) 



r2 



{avM0\S) = 



/•CO i-oo E'^ 

/ dsa{s - ^mD^Ts / dE+e"^^''^ E+ \ - I , (70) 
Jiml J^s V s 
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where the cross section is 

a=^aVl(. + 2m^f . (71) 
3 s-^ 

Numerically solving Eq. (|63p with the above inputs (i.e. numerically solving the 
integrals Eq. (jBH) and Eq. ([701) at each temperature step) , we find thai0 

e 8.5 X 10-1" (^)' , (72) 

where xj \s the final value (T 0) of a:; = T' /T . In Fig. [SI we plot the evolution 
of T' /T, for e = 8.5 x 10""'^°, considering both the numerical solution including the 
effects of the electron mass and the analytic result obtained using Eq. which 
holds in the masslcss electron limit. As expected, the two solutions agree in the 
T < 1 MeV region, where the effects of the electron mass should be negligible. 

In deriving this result we have made several simplifying approximations. The 
most significant of these are the following: a) Using Maxwcllian statistics instead 
of Fcrmi-Dirac statistics to simplify the estimate of {avyi0\£). Using Fermi-Dirac 
statistics should decrease the interaction rate by around 8%. b) We have neglected 
Pauli blocking effects. Including Pauli blocking effects will slightly reduce the in- 
teraction rate since some of the e^' states are filled thereby reducing the available 
phase space. We estimate that the effect of the reduction of the interaction rate 



° For simplicity we have neglected the effect of heating of the photons via e+-e^ annihilations. 
Note that the same effect occurs for the mirror photons which are heated by the annihilations of 
e+'e"', so that Xf is approximately unchanged by this effect. 

0.6 



0.4 



0.3 

i- 

0.2 



0.1 





2 4 6 8 10 

T [MeVJ 

Fig. 6. Evolution of a; = T' /T for e = 8.5 X 10""'^''. The solid line is the numerical solution 
including the effects of the electron mass, while the dashed line is the analytic result obtained 
using Eq. II68I I. which holds in the massless electron limit. As expected the two solutions agree in 
the T > 1 MeV region, where the effects of the electron mass should bo negligible. 
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due to Pauli blocking will be around < 10%. c) We have assumed that negligible 
v'e fi T produced via mirror weak interactions from the e^'. Production of u'^ ^ ^ 
will slightly decrease the T'/T ratio. The effect of this is equivalent to reducing 
the interaction rate by around < 10%. Taking these effects into account, we revise 
Eq. (Ull) to 

e= (1.0 ± 0.10) X 10-9 (^)' . (73) 

Since successful large scale structure studieJ^^l"!^ suggest a rough bound Xf < 0.3, 
then our result, Eq. ([73| . suggests a corresponding rough bound e < 10^^. This 
upper bound is similar to what previously find with the simpler approximation 
of Eq. (|61|) . and again confirm that the current proposed interpretation of direct 
detection experiments in terms of mirror dark matter is compatible with constraints 
coming from early Universe analysis. 



4.3.3. Primordial mirror Helium 



Let us now consider the implications of kinetic mixing for mirror BBN. It has already 
been discussed in previous sections that, compared with the ordinary matter sector, 
we expect a larger mirror helium mass fraction if T' < T, as currently required. 
Essentially, this is because the expansion rate of the Universe is faster at earlier 
times, which implies that the freeze out temperature of mirror weak interactions 
will be higher than that in the ordinary sector. Given our above calculations of the 
T'/T evolution, we can estimate the mirror helium mass fraction as a function of e 
within our model, which we now discuss. 

We follow the same procedure used in Sec. l4.1[ recalling that in the mirror sector 
we have the same relation (|46p than in the ordinary one, except that we change 



t'j^ and Tw 



T^. This is possible since T/v depends only logarithmically 



on the baryon to photon ratio, so that we still have ^ T^r, even though the 
mirror baryon to photon ratio is much larger than the ordinary one. We obtain the 



freeze out temperature from the equality Tw{T{y) 
expressed as follows 



H{T{^), where F and H are 



Tw{T') 
H{T') 



Gl 



(74) 



1.66 g 



1/2 



1/2 



PI 



1.66 



gi/2(l + B-4/3e-8/3r'4/3^ 



1/2 T 



12 



Mpi 



1.66 gi/25-2/3 



-4/3 rp/8/3 



PI 



(75) 



For Eq. ^ we used Eq. ^ with B = 0.74 • 10^^ MeV and considered that 
5-4/3 g-8/3 rpi4/3 jj^m;]^ larger than 1 for e values of interest. Imposing the equality 
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of expressions (|74p and ((75|) we obtain 



^w[T{^) _ G], Afpi 2/3 ,4/3 7./ 7/3 ^ , 
H{T{^) - 1.66 ^ 



(76) 
(77) 



1 /3 

where for the second expression we uscdTw ^ (1.66 g^^^/Gl Mpi) ' , as computed 
from the equation Tw(Tw) = H{Tw)- 

The quantity can be estimated straightforwardly 

^^ -1/2 

(78) 



N 



0.3 g- 



-1/2 



-'at 



where we used In ^ 0.3 g-^^^Mpi/T^ and B-^-^Tn > 1. 

Finahy, using Eqs. (|77p and (|78p. we can estimate the mirror hehum mass frac- 
tion 



2 exp(-t;,/T) 
1 + exp(Am/r; 



2 exp 



-1/2 



1 + exp 



Am 52/7 



.9/7 



.4/7 



w 



(79) 



Using the typical values Tw ^ 0.8 MeV, T/v ~ 0.07 MeV, Am ~ 1.29 MeV, 
B ~ 0.74 • 10^^ MeV, t ~ 886.7 s, In ^ 200 s, we obtain 

2 exp [-2.3 • 10^ e] 2 exp [-2.3 • IQ-^ e^g] ^^^^ 



1 + exp [0.3 • 105 e4/7] ^ 



exp 



0.22 e 



4/7 



Thus, for e ~ 10-^ wc find Y' ~ 0.87. 

In Fig. [7] wc plot the predicted value of Y', together with the ratio T'/T at 1 
MeV, as a function of e, in the interesting region around e ^ 10~^. We see that a 
decrease of the strength of the photon-mirror photon kinetic mixing induces a lower 
mirror temperature, and thus a larger primordial mirror helium abundance (since 
mirror BBN happens at a larger cosmic expansion rate). 

A more precise computation of the primordial mirror helium abundance can 
be done using the calculation of the T'/T evolution in Sec. 14.3.21 and numerically 
solving the equations for the rates of reactions governing the equilibrium between 
protons and neutrons. 

We assume the same initial condition, T' ^ T, as in Sec. 14.3.21 with a thermal 
evolution driven by the photon-mirror photon kinetic mixing, which populates the 
mirror radiation just with e*' and 7'. We may describe the evolution of the ratio 
of the temperatures in the two sectors with the approximate analytical expression 
([68]) . which is vahd for T' > 1 MeV and for T < 100 MeV, where we muhiply 
A by the factor uj w 0.8 in order to take into account the reduction in the rate 
from our approximations (Maxwellian statistics instead of Fermi-Dirac, neglecting 
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Fig. 7. The mirror helium mass fraction versus the photon-mirror photon kinetic mixing strength 
e. The ratio of mirror and ordinary temperatures at T = 1 MeV is also plotted. 



Pauli-blocking factors, massless electrons). Assuming Ti ^ 100 McV, then Eq. 
reduces to 



T 



0.25 



(T/MeV)i/4 Y 10-9 ■ ^^^^ 

In our model the mirror sector starts with a temperature much lower than the 
ordinary one, and later the photon-mirror photon kinetic mixing increases only the 
temperature of mirror electron-positrons and photons, since neutrinos are decou- 
pled. We may thus assume T^i <C T', where T' = Tji ~ Tg', which is a reasonable 
approximation for the e values of interest. In this scenario the only reactions we 
need to consider to compute the He' abundance Y are 



and 



(82) 



We may neglect the neutron decay n' ^ p' + e^' + D' , since the process is much 
slower than the time for primordial mirror nucleosynthesis, that we estimate to 
happen in the first few seconds of the Universe. The reaction rates of the pr ocess es 
(j82p can be adapted from the standard relations present in Weinberg 's book,IIMl m 
which we neglect the Pauli blocking effect on neutrinos because T^/ ^ T' 



An' 



-i-p' 



A(n' + e+' 
\{p' + e-' - 



nOO 

Jo 



n' + iy') = B 



El.pl, dp. 



where 



B 



Glk(l + 3gi)cos^gc 
27r3ft 



+ l]-i , (83) 



(84) 
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Gwk = 1.16637 X 10^^ GcV^^ is the weak coupling constant, measured from the 
rate of the decay process /i+ — > e+ + + v^, = 1.257 is the axial vector 
couphng of beta decay, measured from the rate of neutron decay, and 9c is the 
Cabibbo angle, with cos = 0.9745, measured from the rate of "'^"'O beta decay 
and other 0+ — > 0+ transitions, and Ee'{pe') = {p1' + m^)^/^, E^'{p^') = p,y. For 
n' + e+' p' + v', E^, - Ee' = Q, while for p' + e~' n' + v' , E^' - E^, = Q, 
where Q = ijin — nip = 1.29 McV. The extremals of integrals in Eqs. ([551) arc fixed 
considering that integrations are taken over all allowed positive values of Pe' ■ 

We may rewrite Eqs. ([55)) by substituting q = —E^i for A„'_j.p' and q = E„' for 
Ap'_>.„', thus obtaining 

1/2 

[e-(?+Q)/T' + l]-i^g^ 

1/2 

[e(<7+Q)/T' ^ (85) 
The differential equation for the ratio Xn' of neutrons to nucleons is 

= Ap'_^„'(l — Xn') — Xn'^p'Xn' ■ (86) 

Note that Y' ~ 2Xn' since, as mentioned before, we can neglect n' decay, and thus 
all available mirror neutrons go into He'. 

We have solved the above equations numerically, using the usual time- 
temperature relation for radiation dominated epoch (|57p . where g takes into ac- 
count only the degrees of freedom of ordinary particles, since the contribution of 
mirror particles is negligible given the initial condition T' T. We used the initial 
condition X„/(0) = 0.5 in Eq. ([55[) . and followed the evolution until Xn' reaches the 
asymptotic value. The results are shown in Fig. [51 where we plot the obtained mass 
fraction of mirror helium versus the strength of the photon-mirror photon kinetic 
mixing (e), for the parameter range of interest for cosmology. As expected, we ob- 
tain high values of the primordial He' mass fraction, with Y' > 0.8 for e < 3 x 10~^. 
For the preferred value emerging from the analysis of the DAMA signal, e ~ 10~^, 
we obtain Y' ~ 0.9, which means that the dark matter is largely mirror helium 
dominated. 

4.3.4. Heavier elements 

We can make a rough estimate of the primordial mass fraction of mirror elements of 
carbon mass and heavier. These are produced essentially via three-body interactions, 
the most important of which is the triple alpha process in the mirror sector 



= B 



1 - 



B q\q + Q? 



{q + QY 



^He'+ W+ -^Hc'^ i^c' + y 



(87) 
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Fig. 8. Mass abundance of primordial mirror helium He' versus the strength of the photon- 
mirror-photon kinetic mixing (e). 



The rate for this process can b e ob tained from the rate for the corresponding process 
in the ordinary matter sector and is given by 



•[exp(-0.37T'-i) +30.3exp(-2.4r'~i)] 



where Xc is the mass fraction of C, T' is in McV units, is the mirror baryon 
density in g/cm^ and the rate is in sec^^. He' cannot be produced in significant 
proportion until T' < 1 MeV, which, from Eq. ([81]), imphes that T < 7 MeV 
(assuming T, > 100 MeV). Using pb' w 0.2(T/MeV)3 g/cin, then we estimate 
that the total mass fraction of C that can be produced is expected to be small, 
Xc < 10-8. 

Anyway, this is just the primordial chemical abundance. Of course, light nuclei 
are expected to be processed into heavier nuclei by stellar nucleosynthesis on succes- 
sive populations of mirror stars. In this process of heavy element enrichment of the 
mirror dark matter interstellar medium a crucial role is played by t he high fraction 
of He' inside mirror stars. In fact, accurate studies have showrp2l that the stellar 
evolution is more rapid for a higher initial He' content, and for Y' w 0.9 it can be 
orders of magnitude faster than the standard case of y = 0.25 in ordinary stars of 
the same masses. This means that the enrichment of heavy mirror elements in the 
halo of the galaxy can be plausibly efficient enough to explain the relatively large 
abundances of He', O' and other heavy elements suggested by the direct detection 
experiments 
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5. Evolution of primordial perturbations 

We have shown that mirror baryons could provide a significant contribution to the 
energy density of the Universe and thus they could constitute a relevant component 
of dark matter. Immediate questions arise: how does the mirror baryonic dark matter 
behave? what are its differences from the more familiar dark matter candidates as 
the cold dark matter (CDM), the hot dark matter (HDM), etc.l In this section we 
discuss the problem of the cosmological structure formation in the presence of mirror 
baryons as a dark matter component. We focus on the structure formation theory 
in linear regime, analyzing the trends of relevant scales (sound speed, Jeans length 
and mass. Silk mass) in both sectors and comparing them with the CDM case, and 
finally showing the temporal evolution of perturbations in all the components of a 
Mirror Universe (ordinary and mirror photons and baryons, and possibly CDM). 

5.1. Mirror baryonic structure formation 

We extend the linear structure formation theory (for the standard scenario, see 
Refs . 1 1 26l and [T27|) to the case of dark matter with a no n- negligible mirror baryonic 
component. 

In a Mirror Universe we assume that a mirror sector is present, so that the 
matter is made of ordinary baryons (the only certain component), non-baryonic 
(dark) matter, and mirror baryons. Thus, it is necessary to study the structure 
formation in all these three components. We proceed in this way: first of all, we 
recall the situation for ordinary baryons, and then we compute the same quantities 
for mirror baryons, comparing them with each other and with the CDM case. 

In general, when dealing with the pre-recombination plasma, we distinguish be- 
tween two types of perturbations, namely be twee n "isoentropic" (adiabatic) and 
"entropic" [isocurvature or isothermal) modes ,^1^ while after matter-radiation de- 
coupling perturbations evolve in the same way regardless of their original nature. 

Here we study only adiabatic perturbations, which arc today the preferred per- 
turbation modes, and we leave out the isocurvature modes, whi ch could also have 
a contribution, but certainly cannot be the dominant component .l^^fflSOl^g recall 
that an adiabatic perturbation satisfies the condition for adiabaticity 

Sm = , (89) 

which, defining as usual S = {Sp/ p), relates perturbations in matter (Sm) and radi- 
ation (Sr) components. 

We will now consider cosmological models where baryons, ordinary or mirror, 
are the dominant form of matter. Thus, it is crucial to study the interaction be- 
tween baryonic matter and radiation during the plasma epoch in both sectors, and 
the simplest way of doing it is by looking at models containing only these two 
components. ^^^^ 

According to the Jeans theory,fI2IE22l ^j^g relevant scale for the gravitational 
instabilities is characterized by the Jeans scale (length and mass) , which now needs 
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to be defined in both the ordinary and mirror sectors. Then, we define the ordinary 
and mirror Jeans lengths as 

V t^Pdom V '-'Pdom 

where pdom is the density of the dominant species, and Wg, are the sound speeds. 
The Jeans masses are 

M, = ^Trp, (^^y = ^p, (Aj)^ , M', = inp', , (91) 

where the density is now that of the perturbed component (ordinary or mirror 
baryons) . 



5.1.1. Evolution of the adiabatic sound speed 

Looking at the expressions of the Jeans length ([TO]) , it is clear that the key issues 
are the evolutions of the sound speeds in both sectors, since they determine the 
scales of gravitational instabilities. Using the definition of adiabatic sound speed, 
we obtain for the two sectors 

1/2 /^„f\ 1/2 



dp 
dp 



.1/2 



dp^ 



(w 



/xl/2 



(92) 



where w and w' are relative respectively to the ordinary and mirror equations of 
state p = wp and p' = w' p' . 

First of all, we consider the standard case of a Universe made of only one sector. 
In a mixture of radiation and baryonic matter the total density and pressure are 
respectively p ~ p-y + pt and p ^ p^ = P^i/'^ (remind that pb — 0). Hence, the 
adiabatic sound speed is given by 

1/2 1 / \-l/2 



dp 
dp 



1 

71 



1 



3/Ofc 



(93) 



where we have used the adiabatic condition (|89p. In particular, using the scaling 
laws p,„ cx pojri ( 1 + ^ ) ^ and p^ cx pi^y{l + z)"^ , together with the definition of matter- 
radiation equality ([7]) (where we now consider only baryons and photons), we obtain 



Vs{z) 



1 

71 



l + z, 



cq 



1 



-1/2 



(94) 



In fact, the relation above is valid only for an ordinary Universe, and it is an 
approximation, for small values of x and the mirror baryon density (remember that 
(3 = ri|,/fib), of the more general equation for a Universe made of two sectors of 
baryons and photons, obtained using Eqs. (|93|) and ([8]) and given by 



1 

71 



1 + /? 



cq 



-1/2 



(95) 
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In the most general case, the matter is made not only of ordinary and mirror 
baryons, but also of some other form of dark matter, then the factor 1 + /3 is 
replaced by 1 + /? + /3_da/, where /3_dm = {^m — — /^b- The presence of the 
term [(1 + x*)/(l + /3)] in equation ((95|) is linked to the shift of matter-radiation 
equality epoch (|9]), thus it only balance this effect, without changing the value of 
the sound speed computed using Eq. (|94|) . 

The mirror plasma contains more baryons and less photons than the ordinary 
one, p'lj = ppb and p' = x^pj. Then, using Eqs. ([92|) and ([8]), we have 



V3 V 4.p'J 



V3 



1 



3 / 1 + .T--* \ / 1 + z, 



-1/2 



(96) 



Let us consider for simplicity the case when dark matter of the Universe is entirely 
due to mirror baryons, ilm — ^'b (i-S-, /? ^ 5). Hence, for the redshifts of cosmological 
relevance, z ^ Zcq, we have v'^ ^ 2x^/3, which is always less than Vg ^ l/VS (some 
example: if x = 0.7, v'^ « 0.5 • Vs] if x = 0.3, « 0.1 • Vg). In expression (^5)) it is 
crucial the presence of the factor [(l + a;~'*)/(l + /3~^)], which is always larger than 1 
(given the bounds (jS]) and (j4]) on the mirror parameters) , so that Ug < Vg during all 
the history of the Universe, and only in the limit of very low scale factors, a <C a^q, 
we obtain Ug ~ ~ 1/ Vs. As we will see in the following, this has important 
consequences on structure formation scales. 

Now we define the scale factor corresponding to the redshift 

(1 + zb^) = (abj)-' = ^ = 3.9 • W\nbh^) . (97) 

Since 1 + Zrec — 1100, ordinary baryon-photon equipartition occurs before recom- 
bination only if rtbh^ > 0.026 (which seems unlikely, given its current estimates). 
According to Eq. (|13p . in the mirror sector the scale of baryon-photon equality a'j,^ 
is dependent on x and it transforms as 

"^-= w^^'-^^^"^^- ^^^^ 

If we remember the definition of the quantity Xeq ~ 0.046(J7m^^)~'^, we have that 
for X > Xoq the decoupling occurs after equipartition (as in the ordinary sector for 
fib/i^ > 0.026), while for x < Xgq it occurs before (as for fib/i^ < 0.026). 

Regardless of which sector we are considering, in the radiation era p^ 3> pt,, 
ensuring that Wg ~ 1/ \/3. In the interval between equipartition and decoupling, 
when pb » p-y, Eq. ([93|) gives Vg — y/ ip^/Spb oc a~-^/^. After decoupling there is no 
more pressure equilibrium between baryons and photons, and Vg is just the velocity 
dispersion of a gas of hydrogen and helium, Vg oc a~^. If ilbh^ < 0.026 or x < Xoq 
(according to what sector we consider), photon-baryon equipartition occurs after 
decoupling, and the intermediate situation does not arise. 

It follows that, by taking care to interchange ab-y with a'^^ and Odoc with a'^^^, 
we have for the sound speed the same trends with the scale factor in both sectors. 
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though with the aforementioned differences in the values. The situation whit x > Xgq 
is resumed below 

{const. a < a'fj^ , 

a- <^ < a < a^i,, , (99) 
a-^ a > a'^^^ . 

If we recall that the matter-radiation equality for a single sector (ordinary) 
Universe, (acq)ord, is always bigger than that for a two sectors (mirror) one, (aoq)mir, 
according to 

(1 + x'^) 

(acq)mir = ^ (acq)ord < (acq)ord , (100) 

together with our hypothesis x < 1 (from the BBN bound) and /? > 1 (cosmo- 
logically interesting situation, i.e., significant mirror baryonic contribution to the 
dark matter) , we obtain the useful inequality always verified in a Universe made of 
ordinary and mirror sectors 

a'f,^ < floq < flb^ . (101) 

It's very important to remark (see Ref. I126P that at decoupling drops from 
(p-y/pb) to (pb/pb). Since p-y (x n-yT while pb oc ni,T with (n^/ni,) ~ 10*^ ^ 1 
in the ordinary sector, this is a large drop in Vg and consequently in Aj. More 
precisely, drops from the value {l/3){p^/pb) = (l/3)(^2^/17(,)(l + Zdoc) to the 
value (5/3)(Tdec/™6) = (5/3)(To/mb)(l + Zdcc), with a reduction factor 

Fi{VLbh^ > 0.026) = ^-^^ = 6.63 • lO-^irtbh^) , (102) 

(«s)dcc 

where (Ws)d7c ^^'^ (■'^s)dtc indicate the sound speed respectively just before and 
after the decoupling. If we consider now the mirror sector and the drop in (Wg)^ at 
decoupling, we find for the reduction factor 

F{{x > x^q) = l3x-'^Fi . (103) 

In the case flbh"^ < 0.026, drops directly from (1/3) to (5/3)(Tdcc/mb) = 
(5/3)(ro/m;,)(l + Zdcc) with a suppression 

F2{nbh^ < 0.026) = = 1.9 • 10"^ . (104) 

(«')dcc 

It's easy to find that in the mirror sector the reduction factor in the case a^,^ > a'^^^ 
is the same as in the ordinary one 

F^{X < Xeq) = F2 . (105) 

Some example: for x = 0.7, F{ « 2.9^Fi; for x = 0.5, F{ = for x = 0.3, F{ « 

37/3 Fi. We remark that, if /? > 1, F{ is at least about an order of magnitude larger 
than Fi. In fact, after decoupling (u^)^ = (5/3)(r^^^/mfc) = (5/3)(Tdcc/™b) = (^s)^ 
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(since T'^^^ = Tdoc), and between equipartition and recombination (Wg)^ < (ws)^- 
The relation above means that the drop is smaller in the mirror sector than in the 
ordinary one. Obviously, before equipartition (Wg)^ = (ws)^ = 1/3, and this is the 
reason why the parameter F2 is the same in both sectors. 

In Fig. [9] we plot the trends with scale factor of the mirror sound speed, in 
comparison with the ordinary one. The ordinary model is a typical one with rif,/i^ > 
0.026, while the mirror model has x ~ 0.6 and /? = 2 (this means that mirror 
baryonic density is twice the ordinary one, chosen in these models about four times 
its current estimation in order to better show the general behaviour). In the same 
figure we show also the aforementioned relative positions of the key epochs (photon- 
baryon equipartition and decoupling) for both sectors, together with the matter- 
radiation equality. If we reduce the value Vlbh? , goes toward higher values, while 
Cdec remains fixed, so that for flth^ < 0.026 decoupling happens before equipartition 
and the intermediate regime, where Vs oc a~^^^, disappears. Analogously, if we 
reduce x, a'^^^ shifts to lower values, until for x < Xcq it occurs before the mirror 
equipartition aj,^, so that the intermediate regime for vi. disappears. 

In Fig. [TU] the same ordinary and mirror sound speeds are plotted together with 
the velocity dispersion of a typical non baryonic cold dark matter candidate of mass 
~ 1 GeV. Note that the horizontal scale is expanded by some decade compared to 
Fig.ini because the key epochs for the CDM velocity evolution (the epochs when the 
particles become non relativistic, Onr, and when they decouple, Od) occur at much 
lower scale factors. 




Fig. 9. The trends of the mirror sound speed (ii^) as a function of the scale factor, compared 
with the ordinary sound speed (fs)- The ordinary model has Cjh'^ = 0.08, while the mirror model 
has X = 0.6 and /3 = 2. All the key epochs are also reported: photon-baryon equipartition and 
decoupling in both sectors, and the matter-radiation equality. 
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Fig. 10. The trends of the mirror (d^) and ordinary (vg) sound speed compared with the velocity 
dispersion of a typical non baryonic cold dark matter candidate of mass ~ 1 GeV {vcdm)> ^nr 
and indicate the scale factors at which the dark matter particles become non relativistic or 
decouple. The ordinary and mirror models are the same as in Fig. (9] but the horizontal scale is 
expanded by some decade in order to show the CDM velocity. 



5.1.2. Evolution of the Jeans length and the Jeans mass 

Recalling the definitions ((90|) and (|9ip . and using the results of Sec. l5.1.1] relative to 
the sound speed, we can compute the evolution of the mirror Jeans length and mass 
and compare them with the analogous quantities for ordinary baryons and CDM. 

We find for the evolution of the adiabatic Jeans length and mass of ordinary 
baryons in the case ftbh^ > 0.026 



Aj oc < 



a 
a 

nl/2 



\3 

Mj (X ^ (X { 



a 

const. 

a-3/2 



acq < a < Qbj , 
flb^ < a < fldcc 
Idcc < a • 



(106) 



Otherwise, if ftbh'^ < 0.026, Ob-y > fldcc, and there is no intermediate phase ab-y < 

a < Odec- 

In the mirror sector it's no more sufficient to interchange Ob-y with a[^^ and fldcc 



with 



'dec 



as made for the sound speed, because from relation (jlOip we know that 



in the mirror sector the photon-baryon equipartition happens before the matter- 
radiation equality (due to the fact that we are considering a mirror sector with 
more baryons and less photons than the ordinary one). It follows that, due to the 
shifts of the key epochs, the intervals of scale factor for the various trends are 
different. As usual, there arc two different possibilities, x > Xeq and x < Xcq (which 
correspond roughly to ftbh^ > 0.026 and ilbh^ < 0.026 in an ordinary Universe), 
where, as discussed in Sec. 15.1.11 for the second one the intermediate situation is 



January 13, 2013 14:54 WSPC/INSTRUCTION FILE mir-cosmo-rcvicw 



44 Paolo Ciarcelluti 

absent. 

Using the results of Sec. 15.1.1] for the sound speed, we find the evolution of the 
adiabatic Jeans length and mass in the case x > x^^ 

^2 



Aj oc < 



a 
a 

ai/2 



(X < 
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< a < Ocq 


const. 
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"doc < a • 



(107) 



We plot in Fig. [TT] with the same horizontal scale the trends of the mirror Jeans 
mass compared with those for the ordinary sector; the parameters of both mirror 
and ordinary models are the ones previously used, i.e. flhh^ ~ 0.08, x = 0.6 > Xcq 
and /3 = 2. 

In the ordinary sector the largest value of the Jeans mass is just before decoupling 
(see Ref . I126P , in the interval abj < a < Odcc , where 

Mj(a < fldcc) = 1.47 • 10"A/o (1 + (fifc/i')'" , (108) 

that for a hypothetical ilt — O.lh^^ and /3 = is ~ IO^^Mq. Just after decoupling 
we have 

Mj(a > Odec) = 2.5 • WH'Iq (1 + /3y^/^ {n^h^Y^'^ , (109) 

that for fif, ~ 0.1/i~^ and /? = is ~ lO^il/©. This drop is very sudden and large, 
changing the Jeans mass by Ff^^ ~ 1.7-10""(rJh/i2)3/2. Otherwise, if rj^/i^ < 0.026, 



M 



'J A 




Fig. 11. The trends of the mirror Jeans mass (Mj) as a function of the scale factor, compared 
with the ordinary Jeans mass (Mj). The ordinary model has Q,f,h^ = 0.08, while the mirror model 
has X = 0.6 and /3 = 2. The horizontal scale is the same as in Fig. [9] We remark that the same 
behaviours of the ordinary sector are present in the mirror sector for different intervals of scale 
factor. 
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ab-y > fldcc, there is no intermediate phase a?,^ < a < a^cc, and M]{a < fldcc) is 
larger 

Mj(a < adcc) ^ 3.1 • IO'^^Mq (1 + (nbh^^^^ , (110) 

while after decoupling it takes the value in Eq. (jl09p . so that the drop is larger, 
~ 8.3 • 10""'^^. We observe that, with the assumptions fit = 1 (totally excluded 
by observations) and /3 = 0, A/j,max (which is the first scale to become gravitation- 
ally unstable and collapse soon after decoupling) has the size of a supcrcluster of 
galaxies. 

If we now consider the expression (|98p . we have 

(111) 



1 + /3 



Acq V P / \ 1 + 2 

which can be used to express the value of the mirror Jeans mass in the interval 
Acq < a < a^gj. (where Mj takes the maximum value) in terms of the ordinary 
Jeans mass in the corresponding ordinary interval ab-y < a < adec- We obtain 

M'Aa < a'^^J « (j^^) ^ • Afj(a < adcc) , (112) 

which, for /? > 1 and .t < 1, means that the Jeans mass for the mirror baryons is 
lower than for the ordinary ones over the entire permitted (/3, x) parameter space, 
with implications for the structure formation process. If, e.g., x = 0.6 and /3 = 2, 
then Mj ^ 0.03 Mj. We can also express the same quantity in terms oi ilb, x and 
/?, in the case that all the dark matter is in the form of mirror baryons, as 

< «dec) « 3.2 . 10"A/e r'^'il + p)-^''' (t|^)'^ (^^^')"' ■ (^3) 

If we remember Eq. (|103p , we obtain that for the mirror model the drop in the 
Jeans mass at decoupling is {F[Y^'^ = f3'^/'^x~^/'^{FiY/'^ ^ which, given our bounds 
on X and /3, is larger than {Fi)^/'^. We give here some numerical example: for 
x = 0.7, (F{)3/2 Si 5/33/2 (Fi)3/2; f^^. ^ ^ q g and /3 = 2 (the case of Fig. [n]), 
(Fi')3/2 « 28(Fi)3/2. 

It's important to stress that these quantities are strongly dependent on the values 
of the free parameters x and /3, which shift the key epochs and change their relative 
positions. We can describe some case useful to understand the general behaviour, 
but if we want an accurate solution of a particular model, we must unambiguously 
identify the different regimes and solve in detail the appropriate equations. 

In Fig. [12] we plot the trends of the mirror and ordinary Jeans mass compared 
with those of a typical non baryonic cold dark matter candidate of mass ~ 1 GeV. 
Apart from the usual expansion of the horizontal scale, due to the much lower 
values of the CDM key epochs as compared to the baryonic ones, a comparison of 
the mirror scenario with the cold dark matter one shows that the maximal value 
of the CDM Jeans mass is several orders of magnitude lower than that for mirror 
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baryons. This implies that a very large range of mass scales, which in a mirror 
baryonic scenario oscillate before decoupling, in a cold dark matter scenario would 
grow unperturbed during all the time (for more details see Sec. 15. 2p . 

For the case x < Xeq, both a|,^ and a'^^^ are smaller than the previous case x > 
cCeq, while the matter-radiation equality remains practically the same; as explained 
in Sec. [2l the mirror decoupling (with the related drop in the associated quantities) 
happens before the matter-radiation equality, and the trends of the mirror Jeans 
length and mass are the following 



A'j cx i 



a 

a 



Alj (X — ^ (X < 



a 



1/2 



a 

const. 



a < a'. 



< a < a' 



dec ' 



, < a < Oo 
< a . 



(114) 



In this case we obtain for the highest value of the Jeans mass just before decou- 
pling the expression 

M'Aa < a',^^) « 3.2 • lO^^M^ r'^^l + P)-''^ [j^^''^ i^bh')-' . 

(115) 

In case x = Xeq, the expressions (|113p and (|115l) . respectively valid for x > Xeq 
and X < Xoq, are coincident, as we expect. If we consider the differences between 
the highest mirror Jeans mass for the particular values x — Xcq/2, x = Xcq and 




Fig. 12. The trends of the mirror (Mj) and ordinary (Mj) Jeans mass compared with those 
of a typical non baryonic cold dark matter candidate of mass ~ 1 GeV (Mqdj^i)\ flnr and a,^ 
indicate the scale factors at which the dark matter particles become non relativistic or decouple, 
respectively. The models are the same as in Fig. 1111 but the horizontal scale is expanded by some 
decade to show the CDM Jeans mass, as in Fig. 1101 
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X = 2xoq-, we obtain the following relations 



^J,max(^cq/2) 



15/2 



eg 



3/2 



-^J,max(^oq) 



0.005 A/j^„,,,(xcq) 



(116) 



3/2 



-^J,max(2a;oq) 



1 + (2a;eq)4 



A^j,max(a;cq) « 64 Af^,„,ax(a:cq) . (117) 



In Fig. [13] we plot the mirror Jeans mass for the three different possibilities: x < Xcq, 
X > Xeq and x = Xcq (the transition between the two regimes), keeping constant 
all other parameters. In these three models the matter-radiation equality is the 
only key epoch which remains almost constant. The change in the trends when x 
becomes lower than Xcq, due to the fact that a'^^^ becomes lower than Oeq, generates 
an evident decrease of the Jeans mass. 




Fig. 13. The trends of the mirror Jeans mass for the cases x < Xeq (solid line), x = Xcq (dotted) 
and X > Xoq (dashed). The model with x > Xcq is the same as in Fig. 1111 the others are obtained 
changing only the value of x and keeping constant all other parameters. As clearly shown in 
the figure, the only key epoch which remains almost constant in the three models is the matter- 
radiation equality; the mirror baryon-photon equipartition and decoupling indicated are relative to 
the model with x < Xcq. The change in the trends when x becomes lower than Xeq is also evident, 
due to the fact that a^^^ becomes lower than Oaq. 
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5.1.3. Evolution of the Hubble mass 

The trends of the Hubble length and mass are expressed, as usually, by 

Mh=<<^oc|^,,, (118) 

It should be emphasized that, as for the ordinary baryons, during the period 
of domination of photons (a < a|,^) the mirror baryonic Jeans mass is an order of 
magnitude larger than the Hubble mass. In fact, following definitions we find 

We plot the trends of the Hubble mass in Figs. [15] and [16] together with other 
fundamental mass scales. 



5.1.4. Dissipative effects: collisional damping 

A peculiar feature of the mirror baryonic scenario is that mirror baryons undergo 
the collisional damping as ordinary ones. This dissipative process modify the purely 
gravitational evolution of perturbations. The physical phenomenon is the interac- 
tion between baryons and photons before the recombination, and the consequent 
dissipation due to viscosity and heat conduction. Around the time of recombination 
the perfect fluid approximation breaks down, and the perturbations in the photon- 
baryon plasma suffer from collisional damping. As decoupling is approached, the 
photon mean free path increases and photons can diffuse from the overdense into 
the underdense regions, thereby smoothing out any inhomogeneities in the photon- 
baryon plasma. This effect is known as Silk damping.^l^ 

In order to obtain an estimate of the effect, we follow Ref. 11181 for ordinary 
baryons, and then we extend to mirror baryons. Wc consider the photon mean free 
path 

= ~ lO^^a^X-^ (nbh^y^ cm , (120) 

AeUeCTT 

where is the electron ionization factor, Ue oc a^'^ is the number density of the 
free electrons and ctt is the cross section for Thomson scattering. Clearly, photon 
free streaming should completely damp all baryonic perturbations with wavelengths 
smaller than A-y. Damping, however, occurs on scales much larger than since the 
photons slowly diffuse from the overdense into the underdense regions, dragging 
along the still tightly coupled baryons. Integrating up to decoupling time we obtain 
the total distance traveled by a typical photon 

As - ^ 3.5 {n.h')-'^" Mpc , (121) 

V 5 fldcc 
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and the associated mass scale, the Silk mass, given by 

Afs = l^p, ^ 6.2 X IQi^ {n,hY'^' 



Mr. 



(122) 



which, assuming ftbh^ ~ 0.02, gives Ms ~ 8 x 10^"* M0. This dissipative process 
causes that fluctuations on scales below the Silk mass are completely washed out 
at the time of recombination and no structure can form on these scales. This has 
consequences on large scale structure power spectrum, where small scales have very 
little power. 

In the mirror sector too, obviously, the photon diffusion from the overdense 
to underdense regions induces a dragging of charged particles and washes out the 
perturbations at scales smaller than the mirror Silk scale 

X'sC:i3f{x){f5nbh^)-^/^Mpc , (123) 

where f{x) = x^/^ for x > Xcq and f{x) = {x j XcqY^'^x^l^ for x < Xcq. Thus, the 
density perturbation scales running the linear growth after the matter-radiation 
equality epoch are limited by the length Ag . The smallest perturbations that survive 
the Silk damping will have the mass 

3 



(124) 



which should be less than 2 x 10^^ Mq in view of the BBN bound x < 0.64. 
Interestingly, for x 



Xcq we obtain, for the current estimate of flmh^ and if all the 
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Fig. 14. The trend of the mirror Silk mass (Mg) over a cosmologically interesting range of x, 
which contains Xeq (we considered Qmh^ — 0.15, so Xeq ~ 0.3). The axis are both logarithmic. 
We show for comparison also the values of the ordinary Silk mass (Mg) and of the free streaming 
mass (Mps) for typical HDM and WDM candidates. 
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dark matter is made of mirror baryons, M'g ^ 10^° M©, a typical galaxy mass. 

At this point it is very interesting a comparison between different damping 
scales, collisional (ordinary and mirror baryons) and collisionless (non-baryonic 
dark matter). We know that for hot dark matter (as a neutrino with mass ^10 
eV) Mpg ~ 10^^ M©, for a typical warm dark matter candidate with mass ~1 keV, 
-^FS^*^ ^ — 10^° Mq, while for a cold dark matter candidate with mass ~1 
GeV the free streaming scale is negligibly small, and it has practically no dissipa- 
tion. From Eq. (|124|) it is evident that the dissipative scale for mirror Silk damping 
is analogous to that for WDM free streaming. Consequently, the cutoff effects on 
the corresponding large scale structure power spectra are similar, though with im- 
portant differences due to the presence of oscillatory features, which makes them 
distinguishable one from the other (see next sections). In Fig. [14] we show this 
comparison together with the trend of the mirror Silk mass over a cosmologically 
interesting range of x. 

5.1.5. Scenarios 

After the description of the fundamental scales for structure formation, let us now 
collect all the informations and discuss the mirror scenarios. They are essentially 
two, according to the value of x, which can be higher or lower than Xeq, and are 
shown respectively in Figs. 1151 and [T51 

Typically, adiabatic perturbations for mirror baryons with sizes larger than the 
maximum value of the Jeans mass, which is Afj(aeq) for x > Xe^ and Afj(a^^^) 
for X < Xcq, experience uninterrupted growth. In particular, they grow as Si, cx 
before matter-radiation equality and as (5b cx a after equality. Fluctuations on scales 
in the mass interval Mg < M < Mj.inax grow as Sb oc while they are still 
outside the Hubble radius. After entering the horizon and until recombination these 
modes oscillate like acoustic waves. The amplitude of the oscillation is constant 
before equilibrium but decreases as a~^/^ between equipartition and recombination. 
After decoupling the modes become unstable again and grow as Sb oc a. Finally all 
perturbations on scales smaller than the value of the Silk mass are dissipated by 
photon diffusion. 

Given this general behaviour, the schematic evolution of an adiabatic mode 
with a reference mass scale Afpdt, with Mg < Mpcrt < Mj{acq), is shown in Fig. [TSl 
for X > Xcq. We distinguish between three evolutionary stages, called A, B and 
C, depending on the size of the perturbation and on the cosmological parameters 
flbh^, X and /3, which determine the behaviour of the mass scales, and in particular 
the key moments (time of horizon crossing and decoupling) and the dissipative Silk 
scale. During stage A, i.e. before the horizon crossing (a < Oont < o-cq), the mode 
grows as Sb oa a^; throughout stage B (oont < a < Odcc) perturbation enters the 
horizon, baryons and photons feel each other, and it oscillates; finally, in stage C 
(a > a^cc)' the photons and baryons decouple, and the mode becomes unstable again 
growing as 6b oc a. We remark that fluctuations with sizes larger than Mj(acq) grow 
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uninterruptedly (because after horizon crossing the photon pressure cannot balance 
the gravity), changing the trend from before MRE to a after it, while those with 
sizes smaller than Mg are completely washed out by photon diffusion. 

After decoupling, all surviving perturbations (those with Mpert > -^s) grow 
steadily until their amplitude becomes so large that the linear theory breaks down 
and one needs to employ a different type of analysis. 

If we look, instead, at the schematic evolution of an adiabatic mode with the 
same reference mass scale Mport but for x < Xeq, as reported in Fig. 1161 we immedi- 
ately notice the lower values of the maximum Jeans mass and the Silk mass, which 
are similar. Therefore, for the plotted perturbative scale there are now only the two 
stages A and C. In general, depending on its size, the perturbation mass can be 
higher or lower than the Silk mass (and approximately also higher or lower than 
the maximum Jeans mass), so modes with Mport > Mg grow continuously before 
and after their horizon entry, while modes with A/pert < are completely washed 
out. 

We find that Mj becomes smaller than the Hubble horizon mass Mn starting 
from a redshift!^ 

z'^ ^ 3750 x-^ (nmh^) , (125) 

which is about Zoq for x ~ 0.64, but it sharply increases for smaller values of a;, as 
shown in Fig.[TJ We can recognize this behaviour also watching at the intersections of 




Fig. 15. Typical evolution of a perturbed scale Mpert in adiabatic mirror baryonic dark matter 
scenario with x > Xoq . The figure shows the Jeans mass Mj , the Silk mass Mg and the Hubble mass 
Mh . The time of horizon crossing of the perturbation is indicated by aent . The three evolutionary 
stages are also indicated: during stage A (a < aont < Ccq) the mode grows as S), oc a^; throughout 
stage B (oent < a < a'^^^) the perturbation oscillates; finally, in stage C (a > a'^^^) the mode 
becomes unstable again and grows as Si, oc a. Fluctuations with size smaller than Mg are wiped 
out by photon diffusion. 
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the lines for Mj and Mh in Figs. [T5] and [TBI Thus, density perturbation scales which 
enter horizon at z ~ Zoq have masses larger than M'j and thus undergo uninterrupted 
linear growth immediately after tc^. Smaller scales for which Mj > A/jj would 
instead first oscillate. Therefore, the large scale structure formation is not delayed 
even if the mirror decoupling did not occur yet, i.e. even if a; > Xe^■ 

When compared with non baryonic dark matter scenarios, the main feature of 
the mirror baryonic scenario is that the mirror baryon density fluctuations should 
undergo the strong coUisional damping around the time of mirror recombination, 
which washes out the perturbations at scales smaller than the mirror Silk scale. 
It follows that density perturbation scales which undergo the linear growth after 
the MRE epoch are limited by the length Ag. This could help in avoiding the ex- 
cess of small scales (of few Mpc) in the CDM power spectrum without tilting the 
spectral index. To some extent, the cutoff effect is analogous to the free streaming 
damping in the case of warm dark matter (WDM), but there are important dif- 
ferences. The point is that, alike usual baryons, the mirror baryonic dark matter 
shows ac oustic osci llations with an impact on the large scale structure (LSS) power 
spectrum.l-^Q^t^^tMj j^^ addition, the oscillations of mirror baryons, transmitted via 
gravity to the ordinary baryons, could cause observable anomalies in the CMB an- 
gular power spectrum. This effect can be observed only if the mirror baryon Jeans 
scale Aj is larger than the Silk scale of ordinary baryons, which sets a principal 
cutoff for CMB oscillations. This would require enough large values of x, and, to- 




Fig. 16. Typical evolution of a perturbed scale Mport in adiabatic mirror baryonic dark matter 
scenario with x < Xoq. The value of Mpert is the same as in Fig. 1151 The time of horizon crossing 
of the perturbation is indicated by aent ■ The figure shows the Jeans mass Mj , the Silk mass Mg 
and the Hubble mass A/jj- Unlike the case x > Xeq (shown in the previous figure), now there are 
only the two evolutionary stages A (a < Oent) and C (a > aent). Fluctuations with size smaller 
than Mg are wiped out by photon diffusion, but in this case the Silk mass is near to the maximum 
Jeans mass. 
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gcther with the possible effects on the large scale power spectrum, it can provide a 
direct test for the mirror baryonic dark matter in CMB and LSS observations. (For 
a complete discussion of the CMB and LSS power spectra for a Mirror Universe see 
next sections.) 

Clearly, for small x the mirror matter recombines before the MRE moment, and 
thus it behaves as the CDM as far as the large scale structure is concerned. However, 
there still can be crucial differences at smaller scales which already went non-linear, 
like galaxies. In our scenario, dark matter in galaxies and clusters can contain both 
CDM and mirror components, or can be even constituted entirely by the mirror 
baryons. 

One can question whether the mirror matter distribution in halos can be dif- 
ferent from that of the CDM. Simulations show that the CDM forms triaxial halos 
with a density profile too clumped toward the center, and overproduces the small 
substructures within the halo. Since mirror baryons constitute a kind of coUisional 
dark matter, it may potentially avoid these problems, at least the one related with 
the excess of small substructures. 

Throughout the above discussion, we have assumed that the matter density of 
the Universe is close to unity. If, instead, the matter density is small and a vacuum 
density contribution is present, we have to add that the Universe may become 
curvature dominated starting from some redshift Zcurv Given the current estimate 
f2A — 0.7, this transition has yet occurred and the growth of perturbations has 
stopped around Zcurv, when the expansion became too rapid for it. 

5.2. Evolution of perturbations 

Here we finally consider the temporal evolution of perturbations as function of the 
scale factor a. All the plots are the results of numerical computations obtained using 
a Fortran code originally written for the ordinary Universe and modified to account 
for the mirror sector. 

We used the synchronous gauge and the evolutionary equations presented in 
Ref. 11341 The difference in the use of other gauges is limited to the gauge-dependent 
behaviour of the density fluctuations on scales larger than the horizon. The fluc- 
tuations can appear as growing modes in one coordinate system and as constant 
mode in another, that is exactly what occurs in the synchronous and the conformal 
Newtonian gauges. 

In the figures we plot the evolution of the components of a Mirror Universe, 
namely the cold dark matter^ the ordinary baryons and photons, and the mirror 
baryons and photons, changing some parameter to evaluate their influence. 

First of all, we comment Fig. [TWb'). which is the most useful to recognize the 
general features of the evolution of perturbations. Starting from the smallest scale 



P As non baryonic dark matter we consider only the cold dark matter, which is at present the 
standard choice in cosmology. 
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factor, we see that all three matter components and the two radiative components 
grow with the same trend (as a^), but the radiative ones have a slightly higher 
density contrast (with a constant rate until they are tightly coupled); this is simply 
the consequence of considering adiabatic perturbations, which are linked in their 
matter and radiation components by the adiabatic condition (|89p . This is the sit- 
uation when the perturbation is out of horizon, but, when it crosses the horizon, 
around a ^ 10^^, things drastically change. Baryons and photons, in each sector 
separately, become causally connected, feel each other, and begin to oscillate for the 
competitive effects of gravity and pressure. Meanwhile, the CDM density perturba- 




Fig. 17. Evolution of perturbations for the components of a Mirror Universe: cold dark mat- 
ter (solid line), ordinary baryons and photons (dot-dashed and dotted) and mirror baryons and 
photons (long dashed and dashed). The model is a fiat Universe with Q,m = 0.3, Q.i,h? = 0.02, 
n^h^ = Q 04 (/3 = 2), h = 0.7, X = 0.6, and plotted scales are logA:(Mpc-i) = -0.5 (a) and -1.0 
(h). 
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tion continues to grow uninterruptedly, at first reducing his rate from to In a (due 
to the rapid expansion during the radiation era), and later, as soon as MRE occurs 
(at a 3 X 10^^ for the considered model), increasing proportionally to a. The 
oscillations of baryons and photons continue until their decoupling, which in the 
mirror sector occurs before than in the ordinary one (scaled by the factor x). This 
moment is marked in the plot as the point where the lines for the two components 
move away one from the other. From this point, the photons in both sectors continue 
the oscillations until they are completely damped, while the mirror and ordinary 
baryons rapidly fall into the potential wells created by the cold dark matter and 
start growing as a. We remark that it's important the way in which the oscillation 
reaches the decoupling; if it is expanding, first it slows down (as if it continues to 
oscillate, but disconnected from the photons) , and then it compresses driven by the 




10 10 10 

scale factor (a/a„) 



Fig. 18. The same as in Fig. 1171 but for scales logfc(Mpc ^) = 



-1.5 (a) and -2.0 (6). 
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gravity; if, otherwise, it is compressing, it directly continues the compression and 
we see in the plot that it immediately stops to oscillate. In this figure we have the 
first behaviour in the mirror sector, the second one in the ordinary sector. 

In Figs. [17] and [18] we compare the behaviours of different scales for the same 
model. The scales are given by log fc(Mpc~^) = —0.5, —1.0, —1.5, —2.0, where 
k ~ 2-K / \ is the wave number. The effect of early entering the horizon of small 
scales (those with higher k) is evident. Going toward bigger scales the superhori- 
zon growth continue for a longer time, delaying more and more the beginning of 
acoustic oscillations, until it occurs out of the coordinate box for the bigger plotted 
scale (log/c = —3.0). Starting from the scale logfc = —1.5, the mirror decoupling 




Fig. 19. Evolution of perturbations for the components of a Mirror Universe: eold dark mat- 
ter (solid line), ordinary baryons and photons (dot-dashed and dotted) and mirror baryons and 
photons (long dashed and dashed). The model is a fiat Universe with Vim = 0.3, r2t,/i^ = 0.02, 
Vl'^h'^ = 0.04 (/3 = 2), h = 0.7, x = 0.5 (a) or 0.4 (6), and plotted scale is logfe(Mpc-l) = -1.0. 
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occurs before the horizon entry of the perturbations, and the evolution of the mirror 
baryons density is similar to that of the CDM. The same happens to the ordinary 
baryons too, but for log A: < —2.0 (since they decouple later), while the evolution 
of mirror baryons is yet indistinguishable from that of the CDM. For bigger scales 
(logfc < —2.5), that are not shown, the evolutions of all three matter components 
are identical. 

As previously seen, the decoupling is a crucial point for structure formation, and 
it assumes a fundamental role specially in the mirror sector, where it occurs before 
than in the ordinary one: mirror baryons can start before growing perturbations 
in their density distribution. For this reason it's important to analyze the effect of 
changing the mirror decoupling time, obtained changing the value of x and leaving 
unchanged all other parameters, as it is possible to do using Figs. [TTyb) . [T9l and [20l 




Fig. 20. The same as in Fig. [191 but for x = 0.3 (a) and 0.2 (b). 
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for X = 0.6,0.5,0.4,0.3,0.2 and the same scale logfc = —1.0. It is evident the shift 
of the mirror decoupling toward lower values of a when reducing x, according to 
the law (|lip . which states a direct proportionality between the two. In particular, 
for X < Xcq ~ 0.3 mirror decoupling occurs before the horizon crossing of the 
perturbation, and mirror baryons mimic more and more the CDM, so that for 
X ~ 0.2 the perturbations in the two components are indistinguishable. For the 
ordinary sector apparently there are no changes, but at a more careful inspection 
we observe some difference due to the different amount of rclativistic mirror species 
(proportional to a;^), which slightly shifts the matter-radiation equality. 




scale factor (c/co) 



Fig. 21. Evolution of perturbations for the components of a Mirror Universe: eold dark matter 
(solid line), ordinary baryons and photons (dot-dashed and dotted) and mirror baryons and pho- 
tons (long dashed and dashed). The models are flat with r2m = 0.3, n^h? = 0.02, fij, = AQt (a) 
or (n„ - fli) (no CDM) (b), h = 0.7, x = OA, and logfc(Mpc-i) = -1.0. 
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Obviously, these are cases where the CDM continues to be the dominant form of 
dark matter, and drives the growth of perturbations, given its continuous increase. 
In any case, if the dominant form of dark matter is made of mirror baryons the 
situation is practically the same, as visible comparing Figs. [TW b) and[2lTa') (where 
we see only slight differences on the CDM and mirror baryons behaviours in the 
central region of the plots) , since mirror baryons decouple before than ordinary ones 
and fall into the potential wells of the CDM, reinforcing them. 

Finally, in the interesting case where mirror baryons constitute all the dark 
matter, they drive the evolution of perturbations. In fact, in Fig. [?lT b) we clearly 
see that the density fluctuations start growing in the mirror matter and the visible 
baryons are involved later, after being recombined, when they rewrite the spectrum 
of already developed mirror structures. This is another effect of a mirror decoupling 
occurring earlier than the ordinary one: the mirror matter can drive the growth of 
perturbations in ordinary matter, and provide the rapid growth soon after recom- 
bination, that is necessary to take into account of the evolved structures that we 
see today. 

After all these considerations, it is evident that the case of mirror baryons is very 
interesting for structure formation, because they are coUisional between themselves 
but collisionless for the ordinary sector, or, in other words, they are self-collisional. 
In this situation baryons and photons in the mirror sector are tightly coupled until 
decoupling, and structures cannot grow before this time, but the mirror decoupling 
happens before the ordinary one, thus structures have enough time to grow accord- 
ing to the limits imposed by CMB and LSS (something not possible in a purely 
ordinary baryonic scenario). Another important feature of the mirror dark matter 
scenario is that, if we consider small values of x, the evolution of primordial pertur- 
bations is very similar to the CDM case, but with a fundamental difference: there 
exists a cutoff scale due to the mirror Silk damping, which kills the small scales, 
overcoming the problems of the CDM scenario with the excessive number of small 
satellites. 

6. Cosmic microwave background and large scale structure 

In the last decade the study of the the Cosmic Microwave Background (CMB) and 
the Large Scale Structure (LSS) is providing a great amount of observational data, 
and their continuous improvement provides powerful cosmological instruments that 
could help us to understand the nature of the dark matter of the Universe. In this 
context, the analysis of the CMB and LSS power spectra for a mirror baryonic dark 
matter scenario plays a key role as a cosmological test of the mirror theory. 

In this section we show the CMB and LSS spectra obtained in presence of mirror 
matter, together with the study of their consistence with observational data, that 
provides useful bounds on th e mirror parameter space, stronger than those based 
on the analysis of BBN aloneU^^^ 

The models are numerically computed by means of the same program used for 
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the computation of the temporal evolution of primordial perturbations (presented 
in the previous section), assuming adiabatic scalar primordial perturbations, a flat 
space-time geometry, and different mixtures of ordinary and mirror matter and 
radiation, and cold dark matter. Starting from an ordinary reference model, we 
study the influence of the mirror sector varying the two parameters that describe it 
for a given ordinary sector, replacing a fraction of CDM (or the entire dark matter) 
with mirror matter. Thus, we add to the usual ones two new mirror parameters: 
the ratio of the temperatures in the two sectors x = T' /T and the mirror baryons 
density 17'^, (also expressed via the ratio /? of the baryonic densities in the two 
sectors). The total and vacuum densities are fixed by our choice of a flat geometry: 
f2o = 1 and 17 a = 1 ^ ^m- 

The shapes, heights and locations of peaks and oscillations in the photons and 
matter power spectra are predicted by all cosmological models based on the in- 
flationary scenario, and their features represent specific signatures of the chosen 
cosmological parameters and of the composition of the Universe, which in turn can 
be accurately determined by precise measurements of these patterns. In particu- 
lar, the exact form of the CMB and LSS power spectra is greatly dependent on 
assumptions about the matter content of the Universe. 

A detailed analysis of the dependence and the sensitivity of the mirror CMB 
and LSS power spectra on other cosmological parameters (matter density, ordinary 
baryon density, Hubble parameter, spectral index of scalar fluctuations, number 
of extra-neutrino species) was performed in Refs. [71] and [TH where the author 
compared the dependences on the number of massless neutrino species Ni, and on 
the ratio of temperatures x. As expected, since we remember that the relativistic 
mirror particles can be parametrized in terms of effective number of extra-neutrino 
species, their effects on the CMB are similar, but different for the LSS. 

6.1. The cosmic microwave background for a Mirror Universe 

As anticipated when we studied the structure formation in the previous section, 
we expect that the existence of a mirror sector influences the cosmic microwave 
background radiation observable today; in this section we want to evaluate this 
effect. 

We choose a starting standard model and add a mirror sector simply removing 
cold dark matter and adding mirror baryons. The values of the parameters for this 
reference model are: fio = I, = 0.30, f^A = 0.70, ujb = ftbh'^ = 0.02, n = 1.00, 
h = 0.70, with only cold dark matter (no massive neutrinos) and scalar adiabatic 
perturbations (with spectral index n). This reference model is not the result of a fit, 
but is arbitrarily chosen consistent with the current knowledge of the cosmological 
parameters; however, this is not a shortcoming, because here we want only to put 
in evidence the differences from a representative reference model. 

From this starting point, first of all we substitute all the cold dark matter with 
mirror baryonic dark matter and evaluate the CMB angular power spectrum varying 
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X from 0.3 to 0.7 (around the upper limit set by the BBN bounds). This is shown 
in top panel of Fig. I22[ where mirror models are plotted together with the reference 
model. The first evidence is that the deviation from the standard model is not linear 
in x: it grows more for bigger x and for x < 0.3 the power spectra are practically 
coincident. This is important, because it means that a Universe where all the dark 
matter is made of mirror baryons could be indistinguishable from a CDM model if 
we analyze the CMB only. We see the largest separation from the reference model for 
X = 0.7, but it would increase for hypothetical larger values of x. The height of the 
first acoustic peak grows for x > 0.3, while the position remains nearly constant. For 
the second peak the opposite occurs, i.e. the height remains practically constant, 
while the position shifts toward higher multipoles /; for the third peak, instead. 




Fig. 22. CMB angular power spectrum for different values of x and = compared with a 

standard model (solid line) . Top panel. Mirror models with the same parameters as the ordinary 
one, and with x = 0.3, 0.5, 0.7 and ij^ = Qrnh^ —i^^b CDM) for all models. Bottom panel. Mirror 
models with the same parameters as the ordinary one, and with x = 0.7 and uj'^ = uj),, 2lui,, Swt, iidi,. 
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we have a shift both in height and position (the absolute shifts are similar to the 
ones for the first two peaks, but the height now decreases instead of increasing). 
Observing also other peaks, we recognize a general pattern: except for the first one, 
odd peaks change both height and location, even ones change location only. 

In bottom panel of Fig. [22] we show the intermediate case of a mixture of CDM 
and mirror matter. We consider x = 0.7, a high value which permits us to see 
well the differences obtained changing a;[, from Wb to 4a;b. The dependence on the 
amount of mirror baryons is lower than on the ratio of temperatures x. In fact, the 
position of the first peak is nearly stable for all mirror models (except for a very 
low increase of height for growing oj^), while differences appear for other peaks. In 
the second peak the position is shifted as in the case without CDM independently 
of w^, while the height is inversely proportional to w^, with a separation appreciable 
for < 'ioJb- For the third peak the behaviour is the same as for the case without 
CDM, with a slightly stronger dependence on w^, while for the other peaks there 
is a weaker dependence on w^. A common feature is that the heights of the peaks 
are not linearly dependent on the mirror baryonic density, while their positions are 
practically insensitive to but depend only on x. 

The dependence of the peaks on x and w^, together with other parameters, is 
analyzed in more details in Refs. ITT] and 1741 

6.1.1. The mirror cosmic microwave background radiation 

In the same way as ordinary photons at decoupling from baryons formed the CMB 
we observe today, also mirror photons at their decoupling formed a mirror cosmic 
microwave background radiation, which unfortunately we cannot observe because 
they don't couple with the ordinary baryons of which we are made. Indeed, there is 
in principle the possibility that mirror CMB photons could influence our observa ble 
CMB in case of existence of a tiny photon-mirror photon kinetic mixing,! ^!^ ! -^^ ! but 
its detection would not be possible with present experiments, given the very low 
estimates for the cross section of this interaction. Nevertheless, the study of mirror 
CMB power spectrum is not only speculative, since it is a way to better understand 
the cosmology with mirror matter and our observable CMB. 

We computed four models of mirror CMB, in order to have enough elements to 
compare with the corresponding observable CMBs. The chosen parameter values are 
those usually taken, and the amount of mirror baryons is the same as the ordinary 
ones, while x is taken as 0.9, 0.7, 0.5 or 0.3 in order to explore different scenarios. 
Thus the parameters of the models are: fio = 1, = 0.3, = a;[, = 0.02, x = 0.7 
or 0.5, h = 0.7, n ~ 1.0. In Fig. [23] we plot the ordinary and mirror CMB spectra 
corresponding to the same model of Mirror Universe. 

The first evidence is that, being scaled by the factor x the temperatures in the 
two sectors, also their temperature fluctuations will be scaled by the same amount, 
as evident if we look at the lowest I values (the fluctuations seeds are the same 
for both sectors). Starting from the top- left panel of Fig. [23l we see that the first 
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mirror CMB peak is higher and shifted to higher muhipoles than the ordinary one, 
while other peaks are both lower and at higher £ values, with a shift growing with 
the order of the peak. 

Observing all the panels, we understand the effect of a change of the parameter 
X on the mirror CMB: (i) for lower x values the first peak is higher (for x = 0.5 
it is nearly 1.5 the ordinary one); (ii) the position shifts to much higher multipoles 
(so that with the same horizontal scale we can no more see some peaks). The 
reason is that a change of x corresponds to a change of the mirror decoupling 
time. The mirror photons, which decouple before the ordinary ones, see a smaller 
sound horizon, scaled approximately by the factor x; since the first peak occurs at 
a multipole i (x (sound horizon)"^, we expect it to shift to higher values of ^ by a 
factor x~^, that is exactly what we observe in Fig. [23l 

It's possible to verify that increasing x the mirror CMB is more and more similar 
to the ordinary one, until for x ~ 1 (not shown in figure) the two power spectra 
are perfectly coincident (as expected, since in this case the two sectors have exactly 
the same temperatures, the same particle contents, and then their photon power 
spectra are necessarily the same). 

If we were able to detect both the ordinary and mirror CMB photons, we had 
two snapshots of the Universe at two different epochs, which were a powerful cos- 
mological instrument, but unfortunately this is impossible with only gravitational 
interactions, because mirror photons are by definition completely invisible for us. 
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Fig. 23. Angular power spectra for ordinary (solid line) and mirror (dashed line) CMB photons. 
The models have Oq = 1, = 0.3, oj;, = o;^ = 0.02, h = 0.7, n = 1.0, and x = 0.9 (top-left 
panel), x = 0.7 (top-right panel), x = 0.5 (bottom-left panel), x = 0.3 (bottom-right panel). 
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6.2. The large scale structure for a Mirror Universe 

Given the oscillatory behaviour of the mirror baryons (different from the smooth 
one of cold dark matter), we expect that mirror matter induces specific signatures 
also on the large scale structure power spectrum. 

In order to evaluate this effect, we computed LSS power spectra using the same 
reference and mirror models used in Sec. l6.1l for the CMB analysis. The two panels 
of Fig. [M] show the LSS for the same models as in Fig. [521 In order to remove the 
dependences of units on the Hubble constant, we plot, as usual, on the x-axis the 
wave number in units of h and on the y-axis the power spectrum in units of . The 




Fig. 24. LSS power spectrum in the linear regime for different values of x and u)'^ = Oj^/i^, 
compared with a standard model (solid line). In order to remove the dependences of units on the 
Hubble constant, we plot on the x-axis the wave number in units of h and on the y-axis the power 
spectrum in units of h^. Top panel. Mirror models with the same parameters as the ordinary one, 
and with x = 0.3,0.5,0.7 and uj'^ = flmh^ — (no CDM) for all models. Bottom panel. Mirror 
models with the same parameters as the ordinary one, and with x = 0.7 and uj'f^ = ujf,, 2ui^, SclI;,, 4LJi,. 
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minimum scale (the maximum k) plotted is placed around the limit of the linear 
regime. 

In top panel of the figure we show the dependence on x for different mirror models 
without CDM; in this case, where all the dark matter is made of mirror baryons, the 
oscillatory effect is obviously maximum. The first evidence is the strong dependence 
on X of the beginning of oscillations: it goes to higher scales for higher x, and below 
X ~ 0.3 the power spectrum for a Mirror Universe approaches more and more the 
CDM one. This behaviour is a consequence of the x dependence of the mirror Silk 
scale, that increases for growing x (for details see Sec. 15.1.41 and Refs. l65 | [69H72| : 
this dissipative scale induces a cutoff in the power spectrum, which is damped with 
an oscillatory behaviour (it will be more evident in Figs. [25] and [26l where we extend 
our models to smaller scales within the non linear region). Oscillations begin at the 
same time of the damping, and they are so deep (because there are many mirror 
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Fig. 25. LSS power spectrum beyond the linear regime for different values of x and o;^ = O^h'^, 
compared with a standard model (solid line). The models have the same parameters as in Fig. 1241 
For comparison we also show a standard CDM model with a negligible amount of baryons. 
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baryons) to go outside the coordinate box. In any case the mirror spectra arc always 
below the ordinary one for every value of x. 

The dependence on the amount of mirror baryons is instead shown in the bottom 
panel of the figure, where only a fraction of the dark matter is made of mirror 
baryons, while the rest is CDM. Contrary to the CMB case, the matter power 
spectrum strongly depends on oj^. The oscillations are deeper for increasing mirror 
baryon densities and the spectrum goes more and more away from the pure CDM 
one. We note also that the damping begins always at the same scale, and thus it 
depends only on x and not on oj^, as we know from the expression (|124p of the 
mirror Silk scale obtained in Sec. 15.1.41 

The same considerations are valid for the oscillation minima, which become much 
deeper for higher mirror baryon densities, but shift very slightly to lower scales, so 
that their positions remain practically constant. 

Let us now extend the behaviour of the matter power spectrum to lower scales, 
which already became non linear. Obviously, since our treatment is based on the 
linear theory, it is no longer valid in non linear regime. Nevertheless, even if it 
cannot be used for a comparison with observations, the extension of our models to 
these scales is very useful to understand the behaviour of the power spectrum in a 
mirror baryonic dark matter scenario, in particular concerning the position of the 
cutoff (we recall that its presence could help in avoiding the problem of the CDM 
scenario with the excessive number of small structures). 

Therefore, in Figs. [25] and [26l we extend the power spectra up to k/h = 10 Mpc~^ 
(corresponding to galactic scales), well beyond the limit of the linear regime, given 
approximately hy k/h < 0.4 Mpc^^. 
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Fig. 26. LSS power spectrum beyond the linear regime for two low values of x (0.1 and 0.2) and 
different amounts of mirror baryons (Qj^ = flm — or ^J, = ^CDm), compared witfi a standard 
model (solid line). The other parameters are the same as in Fig. 1241 For comparison wc show also 
a standard CDM model with a negligible amount of baryons. 
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In Fig. [25] we plot in both panels the same models as in Fig.[24l For comparison 
we show also a standard model characterized by a matter density made almost 
completely of CDM, with only a small contamination of baryons {ili, ~ 0.2% instead 
of ~ 4% of other models). In top panel, the x dependence of the mirror power 
spectra is considered: the vertical scale extends to much lower values compared to 
Fig. 1241 and we can clearly see the deep oscillations, but in particular it is evident 
the presence of the previously cited cutoff. For larger values of x oscillations begin 
earlier and cutoff moves to higher scales. Moreover, note that the model with almost 
all CDM has more power than the same standard model with baryons, which in turn 
has more power than all mirror models for any x and for all the scales. In bottom 
panel we show the dependence on the mirror baryon content. It is remarkable that 
all mirror models stop to oscillate at some low scale and then continue with a 
smooth CDM-like trend. This means that, after the cutoff due to mirror baryons, 
the dominant behaviour is the one characteristic of cold dark matter models (due 
to the lack of a cutoff for CDM). Clearly, for higher mirror baryon densities the 
oscillations continue down to smaller scales, but, contrary to the previous case, 
where all the dark matter was mirror baryonic, there will always be a scale below 
which the spectrum is CDM-like. 

An interesting point of the mirror baryonic scenario is his capability to mimic 
a CDM scenario under certain circumstances and for certain measurements. To 
explain this point, in Fig. [^5] we show models with low x values (0.2 or 0.1) and 
all dark matter made of mirror baryons; we see that for x = 0.2 the standard and 
mirror power spectra are already practically coincident in the linear region. If we go 
down to X = 0.1 the coincidence is extended up to k/h ~ 1 Mpc^^. In principle, we 
could still decrease x and lengthen this region of equivalence between the different 
CDM and mirror models, but we have to remember that we are dealing with linear 
models extended to non linear scales, then neglecting all the non linear phenomena 
(such as merging or stellar feedback), that are very different for the CDM and the 
mirror scenarios. In the same plot we also considered a model with x = 0.2 and 
dark matter composed equally by mirror baryons and CDM. This model shows that 
in principle it's possible a tuning of the cutoff effect reducing the amount of mirror 
matter, in order to better reproduce the cutoff needed to explain, for example, the 
low number of small satellites in galaxies. 

These power spectra provide the linear transfer functions, which constitute the 
principal ingredient for the computation of the power spectrum at non linear scales. 

6.3. Comparison with observations 

So far we have studied the behaviour of the photon and matter power spectra 
varying the two mirror parameters, i.e. the ratio of the temperatures of two sectors 
X and the amount of mirror baryons w^. 

Here we compare these models with some experimental data, in order to estimate 
the compatibility of the mirror scenario with observations and possibly reduce the 



January 13, 2013 14:54 WSPC/INSTRUCTION FILE mir-cosmo-rcvicw 



68 Paolo Ciarcelluti 



parameter ranges. 

In the last decades the anisotropics observed in the CMB temperature became 
the most important source of information on the cosmological parameters: a lot 
of experiments (ground-based, balloon and satellite) were dedicated to its mea- 
surement. At the same time, many authors proved that its joint analysis with the 
fluctuations in the matter distribution (they have both the same primordial ori- 
gin) are a powerful instrument to determine the parameters of the Universe. As in 
Sees. 16.11 and 16. 2| we analyze separately the variation of x and oj'^ in the mirror 
models, using now both the CMB and LSS informations at the same time. 

In order to compare our predictions with observations, we report the only avail- 
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Fig. 27. CMB and LSS power spectra for various mirror models with different values of x, com- 
pared with observations and with a standard reference model (solid line) of parameters f2o = 1, 
= 0.25, Ha = 0.75, LUf, = flfjh^ = 0.023, h = 0.73, n = 0.97. The mirror models have the same 
parameters as the standard one, but with x = 0.3, 0.5, 0.7 and oj^^ = fimh^ — i^b (no CDM) for 
all models. Top panel. Comparison of the photon power spectrum with the WMAP and ACBAR 
data. Bottom panel. Comparison of the matter power spectrum with the 2dF binned data. 



January 13, 2013 14:54 WSPC/INSTRUCTION FILE mir-cosmo-rcvicw 



able analysis of t his ki nd, describe d in R efs. [70l [74l where the authors used for the 
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data, and for the LSS t he 2d F survey (in 



particular the binned power spectrum obtained by Tegmark et a The results 

obtained are still valid, even if better data are now available. In order to com- 
pare with the standard CDM results, we use a reference cosmological model with 
scalar adia bati c perturbations and no massive neutrinos with the following set of 
parameteriSSl = 0.25, cut = 0.023, = 0.75, h = 0.73, n = 0.97. As 
usually, we include in this model the mirror sector; for the sake of comparison, in 
all calculations the total amount of matter = ^cdm + f^b + fi'f, is maintained 
constant. Mirror baryons contribution is thus always increased at the expenses of 
diminishing the CDM contribution. 

We start from Fig. I27[ where we assume that the dark matter is entirely due 
to mirror baryons and we consider variations of the x parameter (as in the upper 
panels of Figs. [22] and [24]) . In top panel, we see that with the accuracy of the 
anisotropy measurements the CMB power spectra for mirror models are perfectly 
compatible with data, except for the one with highest x. Indeed, the deviations from 
the standard model are weak for x < 0.5, even in a Universe full of mirror baryons 
(see Sec. 16. 1|) . In lower panel, instead, the situation is very different: oscillations 
due to mirror baryons are too deep to be in agreement with data, and only models 
with low values of x (namely x < 0.3) are acceptable. Thus, we find the first 
strong constraint on the mirror parameter space: models with high mirror sector 
temperatures and all the dark matter made of mirror baryons have to be excluded. 

In Fig. [^Hlwe compare with observations models with the same x, but different 
mirror baryon contents (as in the bottom panels of Figs. and IM)) . The above 
mentioned low sensitivity of the CMB power spectra on doesn't give us indica- 
tions for this parameter (even for high values of x), but the LSS power spectrum 
helps us again, confirming a sensitivity to the mirror parameters larger than the 
CMB one. This is an example of the great advantage of a joint analysis of CMB and 
LSS power spectra, being the following conclusion impossible looking at the CMB 
only. This plot tells us that also high values of x can be compatible with observations 
if we decrease the amount of mirror baryons in the Universe. It is a second useful 
indication: in case of high mirror sector temperatures we have to change the mirror 
baryon density in order to reproduce the oscillations present in the LSS data. 

Therefore, after the comparison with experimental data, we are left with three 
possibilities for the Mirror Universe parameters: 

• high X < 0.5 and low (differences from the CDM in the CMB, and oscillations 
in the LSS with a depth modulated by the baryon density); 

• low X < 0.3 and high (completely equivalent to the CDM for the CMB, and 
few differences for the LSS in the linear region); 

• low X and low a;[, (completely equivalent to the CDM for the CMB, and nearly 
equivalent for the LSS in the linear region and beyond, according to the mirror 
baryon density). 
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Thus, with the current experimental accuracy, we can exclude only models with high 
X and high w^. Our next step will be to consider some interesting mirror models 
and compute their power spectra. 

In Fig. [29] we plot models with equal amounts of ordinary and mirror baryons 
and a large range of temperatures. This is an interesting situation, because the case 
rij, = fib could be favoured in some baryogenesis scenario with some mechanism that 
naturally lead to equal baryon number densities in both visible and hidden sectors. 
These models are even more interesting when we consider both their CMB and LSS 
power spectra. In top panel of figure we see that the temperature anisotropy spectra 
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Fig. 28. CMB and LSS power spectra for various mirror models with different values of mirror 
baryon density, compared with observations and with a standard reference model (solid line) 
of parameters Oq = 1, = 0.25, f2A = 0.75, ujf, = Q.^h'^ = 0.023, h = 0.73, n = 0.97. 

The mirror models have the same parameters as the standard one, but with x = 0.5 and for 
u)'^ = u)fj, 2u)f, , 'iu)fj , iu)f, . Top panel. Comparison of the photon power spectrum with the WM AP 
and ACBAR data. Bottom panel. Comparison of the matter power spectrum with the 2dF binned 
data. 
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are fully compatible with observations until x ~ 0.5, without large deviations from 
the standard case. In bottom panel, we have a similar situation for the matter power 
spectra, with some oscillations and a slightly bigger slope, that could be useful to 
better fit the oscillations present in the data and to solve the discussed problem of 
the desired cutoff at low scales. Let us observe that we are deliberately neglecting 
the biasing problem, given that an indication on its value can come only from a 
fit of the parameters; then, we have indeed a small freedom to vertically shift the 
curves in order to better fit the experimental data. 

Models of a Mirror Universe where the dark matter is composed in equal parts 
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Fig. 29. CMB and LSS power spectra for various mirror models with different values of x and 
equal amounts of ordinary and mirror baryons, compared with observations and with a standard 
reference model (solid line) of parameters Co = 1, Qrn = 0.25, VIa = 0.75, uji, = fli,h^ = 0.023, 
h = 0.73, n = 0.97. The mirror models have the same parameters as the standard one, but with 
Lu'^ = aj(, and x = 0.3, 0.5, 0.6, 0.7. Top panel. Comparison of the photon power spectrum with the 
WMAP and ACBAR data. Bottom panel. Comparison of the matter power spectrum with the 
2dF binned data. 
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by CDM and mirror baryons are plotted in Fig. 1301 Here we concentrate on values of 
X lower than Fig. 1291 because the larger mirror baryonic density would generate too 
many oscillations in the linear region of the matter power spectrum. In top panel 
we show that, apart from little deviations for the model with higher cc, all other 
models are practically the same. In bottom panel, instead, deviations are big, and 
we can still use LSS as a test for models. Indeed, models with x > 0.4 are probably 
excluded by comparison with observations, and this conclusion is valid even if we 
consider the effect of a possible bias that affects the data. Models with lower x are 
all consistent with observations. 




0.01 0.10 

k/h (Mpc') 



Fig. 30. CMB and LSS power spectra for various mirror models witli different values of x and 
equal amounts of CDM and mirror baryons, compared with observations and with a standard 
reference model (solid line) of parameters Co = 1, Clrn = 0.25, CIa = 0.75, uji, = Cli,h^ = 0.023, 
h = 0.73, n = 0.97. The mirror models have the same parameters as the standard one, but with 
ui'^ = uiCDM a-nd X = 0.2, 0.3, 0.4, 0.5. Top panel. Comparison of the photon power spectrum with 
the WMAP and ACBAR data. Bottom panel. Comparison of the matter power spectrum with the 
2dF binned data. 
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Observing the figures, even if it is surely premature (given tlie lack of a detailed 
statistical analysis of mirror models), we are tempted to guess that mirror baryons 
could hopefully better reproduce the oscillations present in the LSS power spectra. 
Thus, we are waiting for a more accurate investigation (in particular for the large 
scale structure) in order to obtain indications for a standard or a Mirror Universe. 

7. Summary 

As in a puzzle representing the whole Universe, researchers working on mirror matter 
are slowly putting together all the pieces, verifying that their positions are the 
right ones, and going on... Proceeding in this way, so far it's emerging a consistent 
picture of the Universe in presence of mirror dark matter, but not all the pieces are 
found yet. The cosmology of mirror matter is more complex than for other usual 
dark matter candidates, since it is self interacting, shows many different physical 
phenomena, and can form a large variety of structures at all scales, similarly to 
what our baryons do. 

The current situation of the astrophysical research in presence of mirror dark 
matter is shown in Fig. 1311 with emphasis on the connections between the different 
theoretical studies and the predicted observable signatures. Solid lines mark what 
is already done (and is compatible with current observational and experimental 
constraints), while dashed ones mark what is still to do. For the last point a special 
importance is covered by star formation and N-body simulation studies, towards 
which people interested in mirror dark matter should address their efforts in the 
near future. 

At present, we have already investigated the early Universe (thermodynamics 
and primordial nucleosynthesis), and the process of structure formation in linear 
regime, that permit to obtain predictions, respectively, on the abundances of pri- 
mordial elements, and on the observed cosmic microwave background (CMB) and 
large scale structure (LSS) power spectra. In addition, we have studied the evo- 
lution of mirror dark stars, which, together with the mirror star formation, are 
necessary ingredients for the study and the numerical simulations of non linear 
structure formation, and of the formation and evolution of galaxies. Furthermore, 
in future studies they should provide predictions on the observed abundances of 
MACHOs and on the background of gravitational waves. Using the results of mir- 
ror Big Bang nucleosynthesis (BBN) and stellar evolution, we have verified that the 
interpretation, in terms of photon-mirror photon kinetic mixing, of the DAMA an- 
nual modulation signal is compatible with the inferred consequences on the physics 
of the early Universe. Ultimately, we will be able to obtain theoretical estimates of 
gravitational lensing, galactic dark matter distribution, and strange astrophysical 
events still unexplained (as for example dark galaxies, bullet galaxy, ...), that can 
be compared with observations. 

I summarize here just the topics described in this review, the reader interested 
in other aspects can refer to the references listed in bibliography. 
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Thermodynamics of 
the early Universe 



Primordial nuclear 
abundances 




Fig. 31. Current status of the astrophysical research with mirror dark matter: soUd lines mark 
what is already done, while dashed ones mark what is still to be done. For an extensive explanation 
see the text. 



A hidden mirror sector with unbroken parity symmetry has no free parameters 
at microscopic level, while macroscopically its presence is parametrized by two addi- 
tional cosmological quantities, the ratio of photon temperatures in the two sectors, 
X = T' /T, and the ratio of baryonic densities, /? = r2'{,/f26. The parameter /? has a 
lower bound /? > 1 in order to be cosmologically relevant, while x < 0.7 is required 
by the current limits on extra neutrino families computed at BBN epoch. 

In presence of a mirror sector that interacts only gravitationally with ordinary 
particles, the evolution of thermodynamical quantities in the early Universe has 
been studied in detail for temperatures below ~ 100 MeV. The equations were 
numerically solved, obtaining the interesting prediction that the effective number 
of extra-neutrino families raises for decreasing temperatures before and after pri- 
mordial nucleosynthesis; this could help solving the apparent discrepancy in this 
number computed at BBN and CMB formation epochs. 

This study is required for primordial nucleosynthesis analyses and numerical 
simulations, suggesting that, if we respect the bound on a;, the presence of the 
mirror particles has a negligible influence on the ordinary BBN. Instead, ordinary 
particles have a big influence on primordial nucleosynthesis in the mirror sector. 
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producing much more helium He' (up to ~ 80%, depending on the value of x) and 
heavier elements (which nevertheless still have a negligible abundanccQ). Such a 
large value of the primordial He' mass fraction will have important consequences 
for the processes of mirror star formation and evolution. 

Considering the interesting possibility that, besides gravity, ordinary matter in- 
teracts with mirror dark matter via photon-mirror photon kinetic mixing of strength 
e ~ 10"^ (this value is compatible with all experimental and astrophysical bound 
applicable to mirror matter), the mirror scenario is able to explain the results of 
current dark matter direct detection experiments. A photon-mirror photon mixing 
of this magnitude is consistent with constraints from ordinary primordial nucleosyn- 
thesis, as well as the more stringent constraints from CMB and LSS considerations. 
In this case the primordial mirror helium He', emerging from BBN in the mirror 
sector of particles, can reach a very large mass fraction, Y' w 90%, with impor- 
tant implications for the mirror dark matter interpretation of the direct detection 
experiments. 

Star formation represents a key ingredient for the galaxy formation, evolution 
and dark matter distribution, but still need to be studied. Stellar evolution, instead, 
has been studied and it was found that stellar lifetimes are 1-2 orders of magnitude 
smaller than ordinary stars with the same masses. This is qualitatively expected, 
since mirror stars have less hydrogen, that is the principal fuel of nuclear fusion 
during most stellar life, while at the same time more helium makes the stars more 
luminous and then more fuel-consuming. A fast stellar evolution means also that 
in the mirror sector the enrichment of heavy elements is faster due to a higher 
supernova rate, and this helps the interpretation in terms of mirror matter of the 
results of direct detection experiments. 

The difference in temperature of the ordinary and mirror sectors plays an impor- 
tant role for the cosmological evolution of the mirror world, since it shifts to earlier 
times the key epochs for structure formation, which proceeds in the mirror sector 
under different conditions. Considering adiabatic scalar primordial perturbations 
as the input, the trends of all the relevant scales for structure formation (Jeans 
length and mass. Silk scale, horizon scale) for both ordinary and mirror sectors, 
were analyzed and compared with the cold dark matter (CDM) case. These scales 
are functions of the fundamental macroscopic parameters x and /?, and in particular 
they are influenced by the differences between the cosmological key epochs in the 
two sectors. The analysis of the temporal evolution of perturbations for different 
values of mirror temperature and baryonic density shows that for x less than a 
typical value Xeq — 0.3, for which the mirror baryon-photon decoupling happens 
before the matter-radiation equality, mirror baryons are equivalent to the CDM for 
the linear structure formation process. Indeed, the smaller the value of x, the closer 
mirror dark matter resembles standard CDM during the linear regime. This study 



The effects of higher primordial mass fractions of heavy elements on the opacity of matter, and 
consequently on the mirror star formation efficiency, still need to be carefully inspected. 
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shows some interesting results: 1) the mirror Jeans mass Mj is always smaller than 
the ordinary one, thus making easier the growth of perturbations; 2) there exists a 
dissipative mirror Silk scale Mg (analogous to the Silk scale for ordinary baryons), 
that for X ~ Xeq has the value of a typical galaxy mass; 3) for x < Xeq we obtain 
Mj ^ Mg, so that all the primordial perturbations with masses larger than the 
mirror Silk masscomel can grow uninterruptedly. 

Differences in the evolution of matter and radiation components translate into 
different signatures in the CMB and LSS power spectra computed for scalar adi- 
abatic perturbations in a flat Universe. Their analysis demonstrated that the LSS 
spectrum is particularly sensitive to the mirror parameters, due to the presence 
of both the oscillatory features of mirror baryons and the collisional mirror Silk 
damping. For x < Xeq the mirror baryon-photon decoupling happens before the 
matter-radiation equality, so that CMB and LSS power spectra in linear regime 
are equivalent for mirror and CDM cases. For higher values of x the LSS spectra 
strongly depend on the amount of mirror baryons. Qualitatively comparing with the 
CMB and LSS data, it was found that for x < 0.3 the entire dark matter could be 
made of mirror baryons, while in the case x > 0.3 the pattern of the LSS power spec- 
trum excludes the possibility of dark matter consisting entirely of mirror baryons, 
but they could be present as a mixture (up to ~ 50%) with the conventional CDM. 
The results are the following: 1) the present LSS data are not compatible with a 
scenario where all the dark matter is made of mirror baryons, unless we consider 
enough small values x < 0.3 ~ Xeq] 2) high values oix, x > 0.5, can be excluded even 
for a relatively small amount of mirror baryons, since we observe relevant effects on 
LSS power spectra down to values of mirror baryon density of the order 17'^ ^ 17b ; 
3) intermediate values of x, 0.3 < a; < 0.5, can be allowed if the mirror matter is 
a subdominant component of dark matter, fi'^ ^ ^ ^cdm] 4) for small values 
of X, X < 0.3, the mirror and the CDM scenarios are indistinguishable as far as 
the CMB and the linear LSS power spectra are concerned. In this case, in fact, the 
mirror Jeans and Silk lengths, which mark region of the spectrum where the effects 
of mirror baryons are visible, decrease to very low values, which undergo non linear 
growth from relatively large redshift. This is an interesting opportunity for mirror 
matter, since for low mirror temperatures we could obtain models completely equiv- 
alent to the CDM scenario at larger scales (when it works well), but with less power 
at smaller scales (when it shows possible problems). Thus, the present situation can 
exclude only models with high x and high fi'f,. 

The results of the linear evolution of perturbations are the inputs for the study 
of non-linear structure formation down to galactic and subgalactic scales, hopefully 
with the aim of N-body simulations, that has still to be done. 

Estimates of the present population of galactic MACHOs and the gravitational 
waves from the mirror sector, that could be observable with the next generation of 
detectors, can be obtained after the study of mirror stellar formation, used together 
with the already existing predistions of mirror stellar evolution. 
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At stellar scales, there is also the possibility that a non-ncgligiblc amount of 
mirror matter is trapped inside stars, and can have observable consequences. They 
have been studied for the peculiar case of neutron stars, for which it's predicted 
that the mass-radius relation is not unique, but depends on the history of each star. 

Concluding, the astrophysical tests so far used show that the mirror scenario is 
viable and consistent at all investigated scales, making mirror particles a promising 
candidate for dark matter. However, we still need to obtain a complete picture of 
the Mirror Universe. 
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